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Abstract: 

The analytic properties of the M = 2 superconformal index are given a physical interpreta- 
tion in terms of certain BPS surface defects, which arise as the IR limit of supersymmetric 
vortices. The residue of the index at a pole in flavor fugacity is interpreted as the index of 
a superconformal field theory without this flavor symmetry, but endowed with an additional 
surface defect. The residue can be efficiently extracted by acting on the index with a dif- 
ference operator of Ruijsenaars-Schneider type. By imposing the associativity constraints of 
S-duality, we are then able to evaluate the index of all generalized quiver theories of type A, 
for generic values of the three superconformal fugacities, with or without surface defects. 
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1. Introduction 

In recent years much has been learned about rigid supersymmetric quantum field theories 
by placing them on curved compact manifolds, notably (squashed) spheres.^ In several cases 

'^More general susy-preserving backgrounds are beginning to be explored [1, 2]. 
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localization techniques have facilitated the exact calculation of the partition functions on these 
backgrounds, possibly in the presence of supersymmetric defects. Four important classes 
of examples are expectation values of circular Wilson loops in Ad J\f = 2 theories on S"^ 
[3], partition functions of 3d A/" = 2 theories on [4, 5, 6], partition functions of A/" = 
1 and M = 2 theories on x 5^ [7, 8] and partition functions of A/" = 2 theories on 
5^ X [9, 10, 11, 12, 13, 14] (the latter are usually referred to as "superconformal indices" 
if the theory is conformal). 

These observables have been indispensable tools in the exploration of the new classes 
of 4(i and ?>d theories, dubbed respectively class S [15, 16] and class TZ [17]^, which arise 
by twisted compactifications of the Qd (2,0) theory on Riemann surfaces and on hyperbolic 
three-manifolds. This line of enquiry has led to the discovery of remarkable Ad/2d and 3d/3d 
relations. For example, according to the AGT conjecture [20], the partition function of a 
class S theory is computed by a Liouville/Toda correlator on the associated Riemann surface, 
whose complex structure moduli correspond the exactly marginal Ad gauge couplings. In a 
similar spirit, the S"^ x partition function of a class S theory, being independent of the 
exactly marginal gauge theory couplings, is computed by a topological quantum field theory 
correlator on the Riemann surface [21]. 

In this paper we present a complete characterization of this 2d TQFT, and thereby 
provide an algorithm to evaluate the superconformal index of the general theory of class S.^ 
The superconformal index is superficially a simpler observable than the S"^ partition function, 
since it can be evaluated by a free field calculation in any theory that admits a Lagrangian 
description. It is still however a very non-trivial function of the fugacities associated to the 
flavor symmetries and of three additional superconformal fugacities {p,q,t). Its expression 
for a generic class S theory, which is non-Lagrangian, is a priori unknown. For a two- 
dimensional slice {q, t) of the three-dimensional superconformal fugacity space an explicit 
description of the 2d TQFT was proposed in [22, 23] and shown to pass many checks. The 
key was to find a complete basis of functions of the flavor fugacity variables that diagonalizes 
the topological structure constants. The diagonalizing functions turn out to be closely related 
to Macdonald polynomials. Here we obtain the general answer, valid for arbitrary {p,q,t), 
by a completely different strategy. The new method sheds light on the appearance of the 
Macdonald polynomials and yields naturally their expected "elliptic" generalization, which 
would have been difficult to obtain in the previous approach. 

■^More precisely, the class TZ consists of A/" = 2 SCFTs which admit an abelian Chern-Simons-Matter 
description, deformed by a certain class of superpotentials. Conjecturally, it includes the 3d theories defined by 
the compactification of the (2,0) theories on 3d manifolds which admit an ideal triangulation. See also [18, 19] 
for related 3d discussions. 

^More accurately we focus on class S theories of type A, that is, the SU{N) generalized quivers. The 
generalization to theories of type-D and E is in progress. 
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The route that will lead us to the derivation of the general index is somewhat indirect, but 
it is guided by some simple physical ideas. A first source of inspiration is the AGT correspon- 
dence, where normalizable Liouville vertex operators are associated with flavor symmetries of 
the Ad gauge theory, while degenerate vertex operators correspond to inserting extra surface 
defects in the Ad theory. The degenerate operators are the key to the solution of Liouville 
theory by the conformal bootstrap [24]: considering their fusion with normalizable vertex 
operators one can derive functional equations that admit a unique solution. By analogy, we 
expect that adding surface defects to the x partition function should correspond to 
inserting special "degenerate" punctures in the 2d TQFT correlator, and that their fusion 
with the ordinary flavor punctures will lead to "topological" bootstrap equations. This is 
indeed what we find. Another useful heuristic principle is that since divergences in a parti- 
tion function must be related to flat bosonic directions, it should be possible to interpret the 
residue of the index at any of its poles in terms of the behavior of the id field theory "far 
away" in moduli space. 

Guided by this intuition, we set out to evaluate the superconformal index of class S 
theories endowed with BPS surface defects.^ This is a very interesting observable in its own 
right. Surface defects are among the least studied objects in four-dimensional quantum field 
theory and only recently they have started to receive proper attention (see e.g. [28, 29, 30]). 
Surface defects are in fact the only defects compatible with localization with minimal {M = 1) 
supersymmetry in four dimensions. While in this paper we focus on M = 2 theories, we hope 
that our results will serve as a stepping stone towards the study of surface defects in M = 1 
theories. 

In theories that admit a Lagrangian description, such as the class S theories of type 

Ai, it should be possible to evaluate the index and the partition function in the presence 

of surface defects by localization techniques. This is a very interesting direction for future 

work. Here we resort instead to a less direct construction that however also applies to the 

non-Lagrangian higher-rank theories. We start with the physical picture of a surface defect as 

the IR end point of a BPS vortex solution. The construction proceeds by embedding a given 

SOFT TiR into a larger theory Tuv such that turning on a spatially constant Higgs branch 

vacuum expectation value (vev) one flows back to the original Tir. If one then modifies the 

UV boundary conditions of the RG flow by considering instead a position-dependent vev, the 

infrared endpoint is Tir endowed with an additional BPS surface defect. The theory Tuv 

has an additional ^^(1)/ flavor symmetry with respect to Tir. We argue that the residues of 

the index of Tuv at some special poles in the U{l)f flavor fugacity capture the index of Tir 

in the presence of surface defects. Although Tuv does not contain any surface defects, the 

''See [25] for a previous attempt to incorporate surface defects in the superconformal index. The index of 
some N = 2 theories in the presence of loop operators has been evaluated in [26] , while [27] studied the index 
in presence of duality domain walls. 
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analytic structure of its index encodes the possibility to generate them by RG flows. By this 
route we are led to formulate a precise prescription to evaluate the index of Tir with extra 
surface defects. 

The prescription can be formulated entirely in terms of the index of Tir- Remarkably, 
adding surface defects to Tir amounts to acting on its index with certain difference oper- 
ators ©(r,s), closely related to the Hamiltonians of the elliptic Ruijsenaars-Schneider (RS) 
model [31, 32, 33]. The difference operators act as shifts of one of the SU{N) flavor fugaci- 
ties. The analogy with AGT is compelling: the action of the difference operator on one of the 
flavor punctures corresponds to the fusion of a degenerate Liouville primary with one of the 
normalizable primaries. Generalized S-duality predicts that one should get the same result 
independently of which flavor puncture the difference operator is acting on. This immediately 
leads to the conclusion that the functions that diagonalize the topological structure constants 
must be eigenfunctions of the difference operators (under the assumption that the spectrum 
is non-degenerate). By this route we are led to a complete determination of the 2d TQFT 
and thus of the index of a general class S theory of type A, with or without surface defects. 
For the two-dimensional slice (g, t) of superconformal fugacity space the eigenfunctions turn 
out to be proportional to Macdonald polynomials and we re-derive the results of [23]. For 
arbitrary values (p, q, t) of the superconformal fugacities the eigenfunctions are not known in 
closed analytic form. As a demonstration of principle that they exist, that they have a non- 
degenerate spectrum and that they are in fact calculable, we discuss a perturbative scheme 
to determine the eigenfunctions for small q or p - an approach that may have independent 
mathematical interest. 

The structure of this paper is as follows. In section 2 we present an RG construction of 
certain BPS surface defects. In section 3 we give a physical interpretation of some special 
poles of the index in terms of such surface defects. We spell out a precise prescription for 
evaluating the index of a theory in the presence of these defects. In section 4 we apply our 
prescription to the case of Ai quivers of class S and recast it in terms of the action of difference 
operators of RS type. In section 5 we study the properties of these difference operators and 
interpret them physically. In section 6 we combine the consistency conditions of generalized 
S-duality and the explicit form of the difference operators to "bootstrap" the index of a a 
general Ai quiver. In section 7 we dimensionally reduce our results to and interpret them 
in the context of 3d gauge theories. In section 8 we present the extensions to the higher-rank 
theories. We close in section 9 with a discussion and a list of open problems. We collect in 
three appendices some additional material not needed on a first reading. In appendix A we 
describe the embedding of the 2d superalgebra that lives at the location of the defect into 
the Ad superalgebra, and recall some results on 2d partition functions. In appendix B we 
present a perturbative approach to the calculation of the elliptic RS wavefunctions. Finally 
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in appendix C we describe the generalization of our results to the case of flavor punctures 
with reduced symmetry. 



2. An RG construction of supersymmetric surface defects 

We begin by discussing a general construction of BPS surface defects, applicable to a large 
class of = 2 superconformal field theories. The construction proceeds by embedding a given 
SCFT TiR into a larger theory Tuv, such that turning on a spatially constant Higgs branch 
vacuum expectation value (vev) one flows back to the original 7/_r. If one then modifies the 
UV boundary conditions of the RG flow by considering instead a position-dependent vev, the 
infrared endpoint is Tir endowed with an additional BPS surface defect. 

In the rest of the paper we will explain how the pole structure of the superconformal index 
can be understood physically in terms of these surface defects. We will need to make some 
assumptions about the flavor charges of the Higgs branch vevs which initiate the RG flow. 
These assumptions, and many of our calculations, can be elegantly stated in terms of gauging 
a U{l)f flavor symmetry of Tuv- There are reasons for which gauging U{l)f should be a bad 
idea: the gauged theory is not UV complete, and has an anomalous U{l)r R-symmetry. We 
believe that these problems are less serious than they appear, and we will sketch how one 
could overcome them. In any case, we also present an alternative analysis that does not rely 
on gauging U{l)f and leads to the same conclusions. 

2.1 RG Row by a constant Higgs branch vev: the gauged perspective 

As our main example, we take Tir to be a generalized superconformal quiver of type Aj\f-i. 
We focus on a link of the quiver, corresponding to an SU{N) gauge group. We cut the 
link and insert an extra node corresponding to a free hypermultiplet in the bifundamental 
representation of SU{N) x SU{N), see figure 1. The resulting superconformal field theory is 
what we call Tuv- Relative to Tir, Tuv has an extra ^^(1)/ flavor symmetry acting on the 
bifundamental hypermultiplet only. We now describe a supersymmetric RG flow that connects 
Tuv and Tir. The flow is initiated by gauging the U{l)f symmetry in the presence of a Fayet- 
Iliopolous (FI) parameter, which introduces a scale and forces the scalars in the hypermultiplet 
associated to the extra node to acquire a vev. We are interested in the symmetry-breaking 
pattern that preserves the diagonal SU (N) of the original SU{N) x SU (N) non-abelian flavor 
symmetry. This is achieved by choosing the vevs to be proportional to the unit matrix, 

Qad = q^aa , Qaa = Q^aa ■ (2.1) 
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TiR Tijv 

Figure 1: Our first RG example. On the left we indicate schematically the quiver for the SCFT Tir, 
and on the right the quiver for Tuv- The dashed lines represent SU{N) gauge groups. 

The non-abelian D-term constraints are then automatically satisfied. The U{1) j moment 
maps are 

/i3 = Tr(|Q|2-|Q|2) /ii + i//2 = TrQQ . (2.2) 

By an 5C/(2)/j rotation (i.e. a choice of an A/" = 1 subalgebra) we align the FI parameters Vi 
along /i3, so that the ^^(1)/ D-term constraints read 

f^3 = N{\q\^-\qf)=v, ^ii + ifi2 = Nqq = 0, (2.3) 

which have a unique solution up to U{l)f gauge rotations. Taking for definiteness f > 0, the 
solution is |gp = v/N, q = 0. 

The U{l)f X SU{N)i x SU{N)2 gauge symmetry is higgsed to a diagonal SU{N), and 
the bifundamental hypermultiplets are eaten up in the process. In other terms, the extra 
node that we added to the generalized quiver is removed by the RG flow and in the IR we 
recover the original theory Tir- A comment about R-symmetry breaking and its restoration 
is in order. The FI parameter breaks explicitly the SU{2)r symmetry of the UV theory to 
an SO{2)ji subgroup.^ The vev further breaks this SO{2)r spontaneously. There is however 
a linear combination 50(2)^ of 50(2)^ and U{l)f which is preserved in the new vacuum. 
In our conventions we assign 

Rq = 1^ fQ = -^^ (2-4) 

so that the linear combination 

R = R+^- (2.5) 

leaves Q invariant. In the IR, one must recover the full SU{2)j^ x U{l)r R-symmetry of the 
N = 2 superconformal algebra. We identify the SO{2)j^ symmetry, which is preserved all 
along the flow, with the Cartan subalgebra of the infrared SU{2)^. 

^However, the theory has still exact Af — 2 supersymmetry, as explained in [34, 35]. 
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Tib Tuv 
Figure 2: Our second RG example. 

We can also apply the same construction to the external leg of the generalized quiver 
that defines Tir- The external leg is associated to an SU{N) flavor symmetry. To obtain 
Tuv-, we break the leg and insert again a free bifundamental hypermultiplet. Again Tuv has 
an extra U{\)f symmetry. One can repeat exactly the same steps as above and define an RG 
flow between Tuv with higgsed U{l)f and Tir- ^ 

While we have focused on two very concrete examples, the basic idea is more general. 
Given an = 2 SCFT with a U{l)f flavor symmetry, we can gauge the U{l)f and turn 
on an FI parameter. The U (1) f symmetry is then higgsed, four real directions in the Higgs 
branch are eaten, and at much lower energies the theory flows to a new N = 2 SCFT Tir. 
In geometric terms, the IR Higgs branch is the hyper kahler quotient of the UV Higgs branch 
by the tri-holomorphic U{l)f isometry which we gauged: we set the moment maps equal to 
the FI parameters, and remove the orbits of the U{1) f action. In order to define the RG flow 
properly, we need to pick a vacuum for the theory, i.e. a choice of vev for the Higgs branch 
fields. With no loss of generality, we can pick a complex structure, or = 1 subalgebra, 
and only turn on a real FI parameter, so that a Cartan subgroup SO{2)r of SU{2)r is not 
explicitly broken. One can then show that there is always a choice of vacuum such that 
S0{2)r is not spontaneously broken, i.e. a point on the Higgs branch where an SO{2)r 
rotation can be compensated by a U{l)f gauge rotation. Indeed, the Higgs branch of the UV 
SCFT is a hyperkahler cone, with a radial coordinate p which is the moment map for 50(2)^. 
We can look for the point of minimum p among solutions of the moment map constraints. 
At that point, the S0{2)r isometry is aligned with the U{l)f isometry, and there is a linear 
combination SO{2)j^ which is unbroken. In general one expects that SO{2)j^ should become 
the Cartan subgroup of the IR SU{2)^ R-symmetry group. 



^There is a small caveat here concerning the difference between SU{N) and PSU{N) flavor symmetry. We 
will come back to it in sec 3.1 



-7- 



2.2 RG flow by a variable Higgs branch vev: the gauged perspective 

We now consider the same setup of figure 1 , but with a more interesting choice of UV boundary 
conditions. We still embed Tir into the larger theory 7c/y, gauge the extra U{l)f flavor 
symmetry and turn on an FI parameter, but rather than choosing the trivial background 
that brings us back to Tir, we consider a position-dependent profile for the fields. A very 
interesting class of backgrounds are the non-abelian vortex solutions that preserve (2,2) 2d 
supersymmetry. Vortices are solutions that depend on two of the four Cartesian coordinates'', 
which we parametrize by a complex coordinate z. 

Vortex solutions localized in the z plane exist in general only if the FI parameter v is 
turned on. Indeed, the tension of a BPS vortex solution is Inrv, where the integer r counts 
the units of magnetic U{1) f flux on the z plane. Away from the location of the vortices the 
scalar fields asymptote to constant vevs selected by the choice of vacuum at infinity. This 
implies that the endpoint of the RG flow is still Tir away from the location of the vortices. 
As we flow to the IR, the tension of the vortices goes to infinity, and they become BPS surface 
defects in Tir. 

The BPS equations for a vortex solution have a simple structure, and can be organized 
into a set of F-term constraints and some D-term constraints. The F-term constraints force 
the complex scalars in the hypermultiplets to be covariantly holomorphic, and to satisfy the 
usual condition fii + ifi2 = for all gauge groups. The D-term equation is deformed to 

^3 = F,,, + V (2.6) 

for all gauge groups, where the FI term v lives in the U{l)f subgroup. 

The first consequence of these equations is that the vev of gauge-invariant chiral operators, 
i.e. chiral gauge-invariant operators in Tijv with zero U{l)f gauge charge, is holomorphic on 
the z plane. As it goes to a constant vev to infinity, the vev is constant everywhere. The 
moduli space of BPS solutions can be identified with the space of solutions to the F-term 
constraints modulo complexified gauge transformations, dropping the D-term equations. One 
can gauge away the anti-holomorphic component of the gauge connection, = 0, so that 
the scalar vevs can be taken to be holomorphic, and actually polynomial in z. One is left 
with a quotient by polynomial complex gauge transformations. 

The moduli space of BPS solutions with a given vortex number r can be rather intricate. 
We wish to focus on backgrounds where the only scalars that receive a vev at infinity belong 
to hypermultiplets of the extra node in figure 1. The baryon operator B = det Q only carries 
U{l)f gauge charge. The only invertible holomorphic U{l)f transformations are constant 

'At this stage we consider gauge theories in Minkowski space. 
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rescalings, so we can write 

r 

B{z) = Pr{z) = l[{z - z,) , (2.7) 

i=l 

up to an overall rescaling. The degree r is the vortex number, and the Zi are gauge-invariant 
parameters identified with the position of the r vortices. 

Although we cannot generally describe the full moduli space of vortex solutions for a 
generic TijVy we can describe a universal subspace, consisting of solutions where only Q 
receives a vev. Then we are effectively considering solutions in M = 2 SQCD with Nj = Nc 
flavors. Supersymmetric vortices in this model have been extensively studied. We refer to 
[36, 37, 38] for reviews and only recall some of the basic facts here. 

In addition to the 2r (real) position moduli, there are 2{N — l)r internal orientational 
moduli, encoded in the holomorphic matrix Q{z) modulo gauge transformations that leave 
detQ(z) = Pr{z) fixed. For example, for r = 1 there is a C X CP^~^ worth of vortex 
solutions: the C factor describes the position, while the CP^~^ factor parametrizes the 
breaking of the diagonal SU{N) gauge group of Tir down to U{\) x SU (N — 1). For general 
r, the vortex moduli space can be described by a U{r) gauged linear sigma model coupled to 
N fundamentals and one adjoint chiral multiplet (which contains the center of mass degrees 
of freedom) [39]. 

As we flow to the IR, we can keep all the r vortices at the origin, and set Priz) = z^. The 
background preserves a diagonal combination j of the angular momentum j in the z plane 
and of the flavor symmetry /, 

J^j + ^f, (2.8) 

which should be identified with the angular momentum of the IR theory. There is also (for 
any choice of Pr{z), in fact) a preserved combination of the R and flavor symmetries, 

R = R+^, (2.9) 

which we identify with the IR R-symmetry. As the vortices have a tension which is propor- 
tional to the FI parameter, they become infinitely heavy surface defects in the IR SCFT Tir. 
We have just described the construction of an infinite family of surface defects 6(^,0) Tir, 
labeled by the number of vortices r. We can repeat the construction by considering s vortices 
in the plane w, orthogonal to the z plane in M^; we consider the background with B{w) = , 
and obtain surface defects in Tir which we label as S(o,s)- Finally we can consider surface 
defects on both planes at the same time, labelled as ©(^ <,). 

We expect that the flat directions corresponding to the transverse positions of the vortices 
(the coefficients {zj} of the polynomial Pr{z)) will flow to r 2d free chiral multiplets in the 
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IR. It should be possible to strip off these free chiral multiplets, and define reduced surface 
defects ©(r,o) ™ Tir. 

While the difference in the holomorphic gauge couplings of the two SU {N) gauge factors 
coupled to the extra node of figure 1 does not affect the bulk IR SCFT, it is a marginal 
deformation parameter of the &[rfi) defects, which behaves like a Kahler parameter for the 
2d (2,2) theory. Thus the 6(r,o) defects should support an exactly marginal twisted chiral 
operator.® 

2.3 RG flow by Higgs branch vevs: the ungauged perspective 

In order to define our RG fiow, we do not really need to gauge U (I)/- We may simply turn on 
the appropriate Higgs branch vev in Tuv ■, and look at physics below the scale set by the vev. 
Let us first consider the constant-vev background. Gauging the U{\) / ate up three scalar 
fluctuations transverse to the moment map level set, and one along the U{l)f orbit. Thus we 
expect to see in the IR the same Tir we encountered in the previous section, together with 
a hypermultiplet which captures these four extra directions. This hypermultiplet transforms 
with unusual quantum numbers under the SO{2)j^ symmetry R = R + //2: it has two 
components of charge 0, and components of charge ±1. The S0{2) ^ is promoted in the IR to 
an SU{2)j^ which acts as the standard R-symmetry of Tir. Under SU (2)^ the hypermultiplet 
scalars transform as a singlet plus a triplet.^ 

In the absence of the U{l)f gauge field, there are no dynamical vortex solutions. But 
we can still consider the RG flow triggered by turning on a position-dependent background 
with B{z) = Pr{z) = in Tuv- Away from z = 0, the theory still flows to Tir times the 
free hypermultiplet. The bulk theory is however modified at the origin by surface defects, 
which we are led to identify with the surface defects defined in the previous subsection. More 
precisely, we expect them to correspond to ©(r,o) (^s opposed to (3(r^o))i since the coefficients 
of Pr are not fiuctuating degrees of freedom in this setup. It is less obvious to describe the 
fate of the free hypermultiplet near the origin. The size of the gauge orbits now scales to zero 
as we approach the origin, as the Higgs branch vev goes to zero there. This suggests that the 
S(^ o) defects are coupled to a surface defect for the free hypermultiplet theory. Consistency 
of this picture with the "gauged" picture of the previous subsection suggests that the degrees 

*We could also construct more intricate variants of the (S(r,o) surface defects, starting from a general linear 
quiver with k SU{N) nodes, and applying our construction to each of the U{1) flavor symmetry groups of the 
bifundamental hypermultiplets. We expect to obtain surface defects with k — 1 exactly marginal twisted chiral 
operators. We leave the exploration of this construction for future work. 

^This is not in contradiction with the standard R-charge assignment for a hypermultiplet, since this hy- 
permultiplet, being free and decoupled from Tir, has an accidental Sp{l) flavor symmetry. The standard 
R-symmetry for the hypermultiplet is recovered as a diagonal combination of S'f/(2)^ and Sp{l). 
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of freedom of the free hyper that survive at the origin should match with the center-of-mass 
motion degrees of freedom of the vortices in the previous description. 



3. The prescription 



In this section, we give a physical interpretation of a class of poles of the superconformal 
index in terms of the surface defects defined in the previous section. We will be led to a 
precise prescription to evaluate the index in the presence of surface defects. 

Let us first review the definition of the superconformal index and fix our basic conventions. 
The index of an = 2 SCFT can be thought of as a trace over states of the theory in the 
radial quantization, i.e. a partition function on x [7],^^ 



X = Tr(-l)^ 



t 

pq 



pjl2 qj34 



Ji 



(3.1) 



We denoted as ji2 as 734 the rotation generators in two orthogonal planes: ji2 = 32 + Ji and 
J34 = 22 — ji with ji^2 being the Cartans of the Lorentz SU{2)i x SU{2)2 isometry of S^. r 
is the U{l)r generator, and R the SU{2)fi generator of R-symmetries. The are fugacities 
for the flavor symmetry generators /j. We will always assume that 



\p\<l 



kl < 1, 



lil < 1, 



1 



pq 



< 1. 



(3.2) 



The particular index (3.1) counts states annihilated by supercharge Q^- (and its Hermitian 
conjugate): this charge has SU{2)ji charge 5, r-charge —5, and SU{2)i x SU{2)2 charges 
(0, —\)- Other choices of the supercharge will give an equivalent index for M = 2 theories [7]. 
Thus the states which contribute to the index (3.1) satisfy [23] 



2 Qi^, Qi^ 



E -2j2-2R + r = 0. 



(3.3) 



Let us also mention here the single letter indices (partition functions) of the basic ingredients 
of A/" = 2 field theories, the hypermultiplet and the vector multiplet. 



Ili^-{p,q,t,a) 



{l-p){l-q) 
p q 



[a + a 



(3.4) 



+ 



pq 
t 



1-p 1-q {l-p){l 



"in this paper we follow the notations of [23]. 
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Here a is the fugacity for the U{1) charge of the half-hypers. To obtain the multi-particle 
index one computes the plethystic exponent the single-letter partition function. For example, 
the full, multi-particle, index of the vector multiplet is 



Xv= PE 



1 



T^\p,q,t) =exp<^^-X^^(p-,g",r 

I _ Tl 



(3.5) 



.n=l 



Given an index I(a, . . . ) of a theory with SU{N) flavor symmetry labeled by fugacities a the 
index of a theory with this flavor symmetry gauged is given by 



N-l 

n 

.1=1 



2TTiai 



A(a) Xy(a) X(a, . . . ) , 



where 



is the familiar SU {N) Haar measure, and 

Xy(a)= PE 
is the index of an SU{N) vector multiplet. 



(3.6) 



(3.7) 



N 



(3.8) 



3.1 Some important poles 

The superconformal index has a rich analytical structure, in particular it exhibits many poles 
in the flavor fugacities. A divergence in the index is associated to the integration over a 
bosonic zero mode in the 5^ x partition function. The introduction of a maximal set 
of superconformal and flavor fugacities is in fact motivated by the desire to regulate these 
divergences. 

We wish to give a physical interpretation of the simple poles of the index and characterize 
their residues. Not to obscure the main point with heavy notations, consider a schematic 
definition of the index, 

I{a,b) = Ti{-lf a-^b^ , (3.9) 
where / and g are some charges. Let us assume that I has a pole in fugacity a. 



The pole is naturally associated to a bosonic operator, O, with charges f = fo and g = go- 
The fact that we have the divergence means that arbitrary high powers of the operator O 
contribute to the index. Although for generic values of the fugacities fluctuations are 
massive and the x partition function is finite, for a^^ = 1 a flat direction opens up 
and the partition function diverges. The residue at the pole is given by I(b~^'^^^'^ , b), so it 
can be written 

i = Tii-lf b<^' . (3.11) 

where g' denotes the shifted charge 

g' = g-^-^f. (3.12) 

Jo 

The shifted charge is preserved in the background where O has a non-zero vev. To characterize 
the residue, we can use the understanding of the divergence as arising from a zero-mode. We 
then expect the divergence of the partition function to be controlled by the behavior of theory 
"at infinity" in the moduli space parametrized by the vev of O, that is, by the properties of 
the IR theory reached at the endpoint of the the RG flow triggered by giving O a vev. We 
should interpret I as the index of the IR fixed of this RG flow. 

By computing the residue of fugacity a we have lost the U{1) symmetry corresponding 
to it in I and thus the IR theory has smaller global symmetry, i.e. giving vev to O explicitly 
breaks some of the symmetries of the UV theory. Note that the number of states contributing 
to I is smaller (or equal) to the number of states contributing to I since after shifting the 
charges some states might cancel out but new states do not start contributing. This is similar 
in spirit to Romelsberger's prescription for computing indices of IR fixed points in the J\f = 1 
setups [8, 40, 41, 1]. 

Let us specialize this general framework to the flavor fugacities poles that arise from 
Higgs-branch zero modes. Consider a general SCFT with a Higgs branch of vacua A4, pa- 
rameterized by a ring of BPS operators {O"}, which sit in spin representations of SU{2)r, 
have charges under the flavor symmetries and carry no angular momenta or r-charge. The 
index receives a contribution from the highest weight component in the SU{2)ji multiplet, 
complex chiral fields which we still denote as 0° . The operator ring matches the ring of holo- 
morphic functions on Ai. For every generator 0°^, all the powers of operator 0°' contribute 
to the index, and summing the geometric series this implies that the index has a pole^^ at 

^'^In the simplest case, this operator generates a ring and all the its powers contribute. More generally, there 
can be relations, which are taken into account by the numerator I{a, b). 

^■^If we specialize to p = g = 0, so that the chiral operators on the Higgs branch are the only operators 
which contribute to the index, the index takes in general the form of a polynomial in t over the product of 
factors {l-t^" O'i'')- In 

some cases, this specialization of the superconformal fugacities relates the index 
to the Hibert series (as defined e.g. in [42]) of the Higgs branch [23]. This relation was explored recently 
in [43, 44, 45]. 
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We now focus on the constant Higgs branch vevs that we have considered in section 2, 
where we have a clear understanding of how the R-charges of the UV theory match the R- 
charge of the effective IR description. We can then hope to relate the residue of UV index at 
the pole with the index of the IR theory. Let us consider the setup of figure 2, and study the 
divergence of the UV index associated to the baryon operator B. The position of the pole can 
be written as t^^a^^ = 1, where a is the fugacity for U{l)f. The 5C/(2)^ Cartan generator 
is R = R — j^f, so if we sit at the pole a = t~^^/ we are simply replacing with t-^, i.e. 
evaluating the index with the fugacities which are appropriate for the IR theory. 

Using Rb = N/2 and = N , the equation t^^^a^ = 1 implies that the index has N 
poles at the locations 

, i = 0,---,N-l. (3.13) 

The residues of all these poles quantify the same physics: the IR fixed point of a theory 
with a vev for baryonic operator. This baryonic operator breaks the U{1) flavor symmetry 
parametrized by o to Zjy symmetry, and looking at different poles of the index corresponds 
to the index of the IR theory with one of the N group elements of Z^v inserted in the trace 
formula. Thus, it makes sense to either sum over the N poles or consider the contribution of 
one of the poles multiplied by N. In our setup of figure 2 the former corresponds to gauging 
the Z^r center of the SU{N) flavor symmetry and thus leaving only PSU{N) as the flavor 
symmetry in the IR: this scenario is the one naturally occuring in the gauged perspective of 
section 2.1. The latter point of view however corresponds to keeping the full SU{N) flavor 
symmetry which is more natural when considering flows initiated by a vev: this will be the 
point of view we will adopt in the rest of the paper. 

In the IR we also expect to have the free hypermultiplet with scalars transforming in the 
1 + 3 irreps of SU{2)j^. The index of such a hypermultiplet captures exactly the divergent 
part of the UV index at t^/^a~^ = 1, as it contains exactly the massless fluctuation which 
causes the divergence. If we start from the index of a standard free hypermultiplet, we can 
get the index of the 1 + 3 hypermultiplet by setting a = t^^^. Up to the divergent zero mode, 
we get precisely Ty^, the inverse of the index of a U{1) vector multiplet. This motivates our 
final prescription: 

I[Tm] =NIv Res„^,i/2 ^ I[Tuv] ■ (3.14) 

The origin of the factor N is explained in the paragraph above. In summary, extracting the 
residue at a = t^/^ in the index of Tuv simply amounts to removing the extra node in the 
quiver (see figure 2) and recovering the original theory 7}_r. This is in perfect agreement with 
the intuition coming from the RG with a constant vev for the baryon operator. 



a = exp 



2m 
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The index also receives contributions from holomorphic derivatives of the baryon oper- 
ator in the 12 and 34 planes, dl2d^iB, so by the same logic as above we expect poles at 
t^Bp^qSg^fn = \ "\Yg wish to relate the residues at these poles to the theories obtained in the 
IR by turning on a space-dependent vev B{z,w) = z^w^. As we have explained in section 2, 
the IR theory is Tir in the presence of an additional surface defect ©(^ <,), with some extra 
decoupled matter coming from the free hypermultiplet. In the 5^ x geometry, the surface 
defects fill the "temporal" and (for r, s 7^ 0) the two maximal circles inside the fixed 
by the ji2 and J34 rotations, respectively (or only one of them for r = or s = 0). 

We are almost ready to compute the index of the theory Tir in the presence of surface 
defects. What remains to be done is to figure out which contribution from the free hypermul- 
tiplet should we strip off. A simple choice is to just strip off the bulk contribution. We will 
later identify that with the index for the theory in the presence of the ©(r,s) surface defect, 




(3.15) 



Again, when we set a = t^/'^p^/^ q^l^ we are shifting the UV rotation generators to become 
the rotation generators appropriate to the IR SCFT, as explained in section 2. 



3.2 A toy model: the free hypermultiplet 

To understand which terms should be stripped off from the index to get rid of the decoupled 
degrees of freedom, it is useful to look at a toy-model where the UV theory coincides with a 
single hypermultiplet. The index of a free hypermultiplet is 

1 - t-l/2pn+l^m+l^±l 
1 - t- 

n,m>0 



n .r(t.a";,.p). (3.16) 



Here and throughout the paper we use a shorthand notation where the it exponent means 
that we take the product over both choices of sign. This index is thus given by a product of 
two elliptic Gamma functions. The denominator of (3.16) comes from the scalar fields in 
the hypermultiplet and their derivatives. If we specialize to a = t^/'^p^q^ in the free hyper- 
multiplet index, remove the zero mode and multiply it by Xy, we get some very suggestive 
results. For example, at a = t^^'^p 



^ yr {l-t-^q^+^){l-tq-^) ^ e{t;q) 
""^'^ li^ (1 -p-i^-)(l -pg-+i) eip-^;q) ' '> 



^^The relevance of these functions to the index computations was observed in [41]. See also [46] for a nice 
review of the properties of elliptic Gamma functions. 
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Here the theta-function is defined as 

oo 

e{z- g) ^ J] (1 - zq^){l - z-\^+^), e{qz; q) = 9{z-';q) . (3.18) 

m=0 

We will also use the Pochhammer symbol, 

oo 

(a; h) = JJ(1 - a 6^) . (3.19) 

The residue Ri^ takes the form of a 2d index for a free 2d chiral multiplet, living on the 
circle fixed by the ji2 rotations. We present the details of the precise relation between (3.17) 
and the 2d index in appendix A. The bosonic scalar field in the chiral multiplet has charge 
1 under rotations in the 12 plane: it is exactly the 2d chiral field which we planned to strip 
off from ©(1,0) to get ©(i,o)- 

In order to read off the residue at a = t^/^p^, we can simplify some manipulations by 
looking at the single- letter partition function. We specialize to a = t^/'^p^ , subtract 1 to 
remove the pole 

p q , pq/t-t ^ ip" - p'^^^q + p~'' - 1~ ^ ^^^q) 



1-p 1-q {l-p){l-q) {l-p){l-q) 
or 

l-p'pq-t ^ I - p^"" p~^ - qr^ 
1—p l — q 1 — p~^ 1 — <? 

which gives the residue 

^ 1 / 1 I — 1 — II t'n-\-^ \ / 1 I ')! m \ 1 



This has the form of a 2d index for r free chiral multiplets, living on the circle fixed by the 
ji2 rotations (see appendix A for the definition of the 2d index). The 2d chiral fields we see 
in Rrfi carry charges 1, • • - r under the 12 rotations, and so look exactly like the 2d chirals 
which we planned to strip off from &(r,o) to S^t ©(^.^o)- 

Finally, in order to read off the residue at a = t^^'^p^q^ we can use again the single-letter 
partition function 

P q , Pq/t - t ^ tp^'q' - + p-''q-' - t'^p^-^'q^-' _ ^ 



1-p 1-q {l-p){l-q) (l_p)(l_g) 

We can write that as a sum of three terms. The first term corresponds to i?r,o- The second 
to Rq^s, which is the same as Rsfi with p and q exchanged, 

^ n IT ~ - ^g^"') ^ n ^(^^"'^^ f3 94^ 
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The third piece is 

(t-P5)^f^'f^-r'-p-V-')i^^ (3.25) 

1 — p 1 — q I— p^l— 

This contributes to the residue through rs ratios of the form 

(l-p"+lg'"+l)(l-t-V"g-"-) ^pq 

{I - p-n-lq-m-l-^^l _ ^pUgm-^ i ■ > 

Hence the total residue is 

^r,s = (^) RrfiRo,s (3.27) 

Up to the prefactor, this has the form of the product of two 2d indices, one for a 2d theory on 
the circle fixed by ji2, the other living on the circle fixed by 734. The prefactor corresponds 
to a re-definition of the U{l)r charge by a constant shift. Such shifts of an abelian symmetry 
generator are common in the presence of defects, unless the generator sits in a non- abelian 
factor of the subgroup of the superconformal group preserved by the defect. Here we have 
two distinct surface defects, in a quarter-BPS configuration. 



Inspired by this calculation, we propose our final recipe to completely strip off the free 
hypermultiplet contributions: 



N R I Iv Res^=(i/2.„r, 



1 



pr/N qa/N 



uv\ 



(3.28) 



3.3 Gauging U{l)f and the index 

Our basic prescription (3.28) has a suggestive interpretation in terms of the index of the 
theory where the U{l)f symmetry is gauged, 

f^a^Iv^iTuv]. (3.29) 

Here the contour integral is on the unit circle |a| = 1 and we added a FI parameter ^ (in units 
where the has unit radius), which can be turned on in the index. 

Notice that (3.29) only makes sense if ^ is quantized. If > 0, it is natural to try and 
evaluate the contour integral by picking the residues at |a| < 1. In particular, we will pick 
the residues at a = t^^'^p''^^ q^/^ . In general there could be other poles. The dominant 
contribution at large ^ will be the residue at a = t^^^, which is precisely our prescription for 

^*We thank N. Seiberg for comments on this point. 
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the index of 7}_r, while the subleading contributions (or a subset of them) correspond to the 
index of 7/_r with the insertions of &(^r,s)- The UV index can then be written as a sum over 
the indices of the IR theories associated to the different supersymmetric backgrounds. 

One may worry that the U{l)f gauged theory is not well defined since it has both a 
Landau pole and an anomalous U {l)r symmetry. The presence of a Landau pole may not be 
a serious problem, as the index is independent of the gauge coupling, and we can suppress 
the Landau pole arbitrarily by making the gauge coupling smaller and smaller. The anomaly 
in U{l)r would appear to be a more serious obstruction. In principle, one can reabsorb the 
TtF a F anomaly for the U{l)r symmetry by redefining the U{l)r charge as 

r^r + bo AaF, (3.30) 

where 60 is some proportionality constant. The Chern Simons action is usually not gauge 
invariant, but for an Abelian gauge theory (F = dA) C/(l)j, on a three manifold with no 
two-cycles, such as S"^, the gauge variation 

/ dAAF (3.31) 
Js3 

is actually zero. 

The redefinition of the r charge will place an overall power oi pq/t in front of the con- 
tribution to the index of interesting topological sectors for U{l)f. A pair of linked vortices 
may exactly fit the bill. Each vortex contains a unit of flux in the core, i.e. the integral of 
A around the vortex should be one. If we have linked vortices, at the core of each we have 
flux F from that vortex, and potential A from the other vortex, which conspire to give a unit 
contribution to J^a A A F. This may neatly explain the prefactor {^Y^ we encountered in 

Rr,s- 



4. Residues and Ruijsenaars-Schneider models: Ai theories 

We now turn to a detailed analysis of the analytic properties of the index for the Ai theories of 
class 5, the Ai generalized quivers of [15]. The generalization to higher rank is conceptually 
straightforward and will be presented in section 8. The Ai case is technically the simplest 
as there is only one type of flavor puncture, carrying an SU (2) flavor symmetry. The basic 
building block is the free half-hypermultiplet in the trifundamental representation of SU{2)^ , 
associated to the three-punctured sphere. The general Ai theory of class S is associated to a 
pair-of-pants decomposition of a genus g surface with p punctures. 
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Figure 3: The setup of figure 2 for the Ai generahzed quivers. The SCFT Tir, associated to 
a Riemann surface C with s punctures, is coupled to the SCFT of free trifundamental half-hyper, 
associated to a sphere with three punctures. The two SCFTs are coupled by choosing an SU{2) flavor 
symmetry of each theory, and weakly gauging a diagonal SU{2) subgroup. The resulting SCFT is 
denoted by 7c/y, and corresponds to the degeneration limit of a surface with s + 1 punctures shown 
in the figure. 

We consider the setup of figure 2. The SCFT Tuv is taken to be an Ai class S theory 
with p = s + 1 punctures, s > 1. We are interested in the degeneration limit where a three- 
punctured sphere is connected by a long cylinder to the rest of the surface, an s-punctured 
Riemann surface which we denote by C, see figure 3. The field theory interpretation of this 
geometry is familiar: the SCFT associated to C is coupled to a trifundamental half-hyper, by 
weakly gauging a diagonal SU{2) group. In the notation of figure 2, the theory associated to 
C is Tir. 

We then focus on the functional dependence of the superconformal index of Tuv either 
one of the flavor fugacities of the weakly-coupled hyper. As outlined in sections 3 and 4, we 
expect to find poles at certain values of the flavor fugacity, and to be able to interpret the 
residue at each of these poles as the index of Tir in the presence of additional surface defects. 
The analysis will be somewhat technical but the main idea is simple. The basic point is that 
since the index of Tuv is a contour integral in the fugacity of the gauged SU{2) associated 
to the long cylinder, a pole in the flavor fugacity arises when the integration contour is 
"pinched" by two singularities. A little calculation reveals that the residue takes the form 
of a difference operator acting on the index of Tir- Our main claim is that the action of 
this difference operator introduces a surface defect in Tir, equation (4.18). Remarkably, the 
difference operator turns out to be the Hamiltonian of the relativistic Calogero-Moser model, 
also known as the elliptic Ruijsenaars-Schneider (RS) model. 

^^The RS integrable models appear in several contexts in the study of A/" = 2 gauge theories, see e.g. [47, 48, 
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In the following sections we will employ these difference operators to extract several 
physical statements about the index of theories of class S. In particular we will gather 
evidence for the interpretation of the difference operators as operators introducing surface 
defects into the index computation. The difference operator can act on any one of the s 
flavor fugacities of Tir, and by generalized S-duality we expect to obtain the same result 
regardless of which fugacity we decide to act on. This constraint is so powerful that it allows 
to completely fix the index for any Ai generalized quiver, with or without surface defects. 
Of course, in the Ai case all theories have an explicit Lagrangian description, and the index 
of any quiver admits an explicit representation as a matrix integral. What this new logic 
buys us is a representation of the index where the structure constants of the associated 2d 
TQFT have been diagonalized. More importantly, this algorithm will generalize easily to the 
higher-rank quivers, which generically do not admit a direct Lagrangian description. 

For general values of the three superconformal fugacities {p,q,t), the answer is given 
in terms of eigenfunctions and eigenvalues of the elliptic RS model, whose analytic form 
is not fully known. We construct however approximate eigenfunctions in a certain limit, 
which serves as a proof of concept of their existence and uniqueness (see appendix B). For 
p — )■ (or equivalently g — )• 0), the elliptic RS model degenerates to the much better studied 
trigonometric RS model, whose eigenfunctions are the Macdonald's polynomials. By this 
route we recover from "first principles" the results of [23] . Finally we consider the dimensional 
reduction of the 4d superconformal index to the 3d partition function. 



Let us introduce the basic ingredients. The index associated to the three-punctured 
sphere, the index of the free trifundamental half-hypermultiplet, is given by 

Ihyp{a,b,c)= n ^ ^ r— — = r(t-2a^^b^^c^';q,p) . (4.1) 

m,n>0 ^ p q l^a C 

We couple this free SOFT to a general Ai quiver SCFT, which we call Tib, associated to 
a Riemann surface C, by gauging fugacity c. The result is the SCFT that we call 7[/y, see 
figures 2 and 3. The index I[Tuv] (= ^ for short) reads 

/dc 
—Aic)Iv{c)Ihypia,b,c)Ic{c-\ ...), (4.2) 

49, 50]. The same quantum mechanical integrable models are also related to gauge theories in lower dimensions 
(see e.g. [51, 52, 53, 54, 55, 56]). It would be interesting to understand the connection between these lower- 
dimensional gauge theories and the 2d TQFT that computes the index, which for q = t is known to reduce to 
2d q-deformed Yang-Mills in the zero area limit [22]. 
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where Tq = I[Tir]- From the physical considerations of sections 2 and 3, we expect X(a, . . . 
to have simple poles at a = t^/^p'^/^g*/^, for r and s non-negative integers. 
Let us start with the simplest case, the pole at 



1 



a = t2. (4.3) 

To verify that I has indeed a pole at this value of the fugacity, we proceed as follows. The 
integrand has (among others) simple poles in c at 

Setting a = t2 these two poles collide and pinch the c integration contour producing a simple 
pole in a. To compute the residue at o = of X we can pick up the residues at the poles 
(4.4) in c inside the integration contour. The other poles in c do not contribute to the residue 
in a. Explicit evaluation gives 

Res 1 - [XI =y^Res i i , - \A(c)Iv(c)Ihvv(0',b,c)Ic(c~^,...)] 



1 l2Zc{h,...)A{b)Zv{b) 



PE 



(4.5) 



2 2 n±(i-^^') 

The origin of the different numeric factors is as follows. The factor of 2 comes from summing 
over the two poles in (4.4). One factor of ^ comes from the Haar measure of SU{2) and the 
second factor of ^ comes from evaluating the hypermultiplet index at the pole. 

There is also a pole at a = — ta, whose residue is \Xy^ Tc{—h, . . .). As we discussed in 
section 3 we need only consider the pole (4.3) since with either sign the poles describe the 
same physics. 

All in all, we can summarize our calculation by the simple relation^^ 

1[Tir] = 2 Tv Res _ i - [T[Tuv]] , (4.6) 



^®We do not claim these are the only poles in general. The poles we discuss here are the ones for which we 
found a simple interpretation in terms of the behavior at infinity of the Ifiggs branch. Other poles may have 
different interpretations; see e.g. [43] for some special examples. 

^''Let us reiterate in this concrete context the logic of section 3.f : in the gauged perspective, as the contour 
integral corresponding to the gauging of f/(l)/, is around the unit circle both poles at ±45 contribute. The 
index of the IR theory is invariant under negating any of the flavor fugacities. In other words the flavor 
symmetry has additional discrete Z2 component which is gauged when we compute the residues. Thus, after 
gauging the U{l)f flavor symmetry (which for Ai quivers is enhanced to 5(7(2)), we only have PSU{2) flavor 
symmetry, as opposed to SU{2), associated with one of the punctures. In the un-gauged perspective it is more 
natural to consider one of the poles, say at a = , taken twice and thus keeping the center of 5(7(2) not 
gauged. 

^^Recall the symbols I[Tir] = Ic are used interchangeably. 
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which precisely confirms (3.28) for r = s = 0. Taking the residue of the UV index at the 
"extra" U{l)f flavor fugacity a = we have "completely closed" the U{l)f puncture and 
got back the IR index. 



Let us now consider the poles with non-trivial p and q dependence, 

a = ±t'2 f/^ q'/\ (4.7) 

with non-negative integers r and s. We need only focus on the positive sign. We give the 
details of the calculation for s > 0, r = 0, as the generalization to non-zero r is straightforward 
we will simply quote the result at the end of this subsection. We start by considering the 
poles of the integrand of (4.2) at 

c = tlg'"ii-, c-^=t^2q^2-^ (4 8) 

a 



a h ' 
h 

o ' ^ ah 

If the non-negative integers rrii are chosen such that 

s = mi + 771,2 1 (4-9) 



1 1 , 



then the pairs of poles in each line of (4.8) collide when a = ga*. Further, since all 
expressions are symmetric under c — )• the residues at the poles corresponding to the two 
lines in (4.8) are equal under exchange of mi with m2. It follows that to compute the residue 
at (4.7), we need to keep track of the terms in the c contour integral classified by a partition 
of s into two non- negative parts. E.g. for s = 1 we get two different choices, mi = 1, m2 = 
and mi = 0, m2 = 1. Let us then compute the residues. The contribution of the hyper is 
given by 



^w=2^xc(<zi--6,) n n ' ' 



i,j=lm,n>0 ^ ^ ^ ^ J' ^ 



(4.10) 



where we have defined hi = 6, 62 = and the prime over the second product indicates that 
we are omitting the diverging term. The factor of 2 in front is coming from the two lines 
in (4.8), while the factor of half arises since the pole in a appears as 1/(1 — tq^a~'^). 

We have to multiply each of these factors by the index of the vector multiplet and the 
Haar measure evaluated at the pole and sum over all the partitions of s, which we will denote 
by 7r(s). The index of the vector multiplet together with the Haar measure evaluated at the 
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pole is given by 



^1 

ij=l m,n>0 

(4.11) 



Combining the two factors (4.10) and (4.11), and further multiplying by twice the index of 
free vector multiplet, we finally have 



(4.12) 



{mi}e7r(r) {mi}G7r(r) 



We see that the computation of the residue at the pole (4.7) amounts to applying a difference 
operator S(o,s) to the index, 

©(0,,)Xc^2XyRes 1 ^ [X] . (4.13) 

The generalizations of this calculation to non-zero r in (4.7) is immediate. The residue 
is given by a sum over terms coming from partitions of both s and r, 

2 2 

s = ^^mi, r = ^^m^. (4-14) 

i=l 1=1 

We will denote the operator which computes the residue at (4.7) by &(^r,s)- 



by 



(s=l) 

Let us now compute explicitly the basic operator ©(o,i)- First the factors /j^.j are given 

(,= 1) ^ yr {1 - p"^t){l - p"^+H~^) {l-p^+\t~^h2/hl){l-p^q~Hbl/h2) 

(1 -p'"+i9)(l -p"*?-!) (l-p'^+^hi/h2){l-p'^h2/hi) 
9{t-p) 0C-b,/b2;p) ^^^^^^ 



0{q-';p) 0{b2/bi-p) • 
Thus the residue at a = t^q2 is given by the action of the following operator 

2—1 



9{t;p) 



-1. 



;p) 



e{b^-p) '■■■'^ 9ib~^;p) ' 
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The operator 6{o,i) is, up to a conjugation by a simple function (see appendix B), the basic 
difference operator of the Ai RS model! We recognize the prefactor 



e{t;p) 



as the function 



i?o,i introduced in (3.24). This is the 2d index of one free chiral field in two dimensions (see 
appendix A) and we interpret it as the center of mass degree of freedom of the surface defect. 
As we have discussed in section 3, it is natural to strip off this factor from the definition 
of the defect. We can then define the difference operator associated to the "bare" defect, 
6(0,1) = -^0 i®(o,i)- However, we will mostly write equations for the "full" operators S(r,s), 
since they are somewhat easier to manipulate algebraically. 

The operator corresponding to residue at a = t2p2^ ®(i,o)) is simply obtained by ex- 
changing q and p in (4.16). It is not difficult to write down the operators for general s, 



e{tq^-p) 



m—l 

n 

.n=0 



9{tq 



s~2m+nu2 . 



b';p) 



—m—l 

n 

n=0 



-^■p) 



-m+nl)2.p-^ 



Finally we quote the general result for all r and s, 



(4.17) 



^[TlR, &{r,s)] = 2 Xy Res 1 



I [I[Tuv]] = 6(r,s) ATir] 



Ei=i"i='-Ei=i' 



n 



'■nii-l' 

n 

m=0 



e{p-^iq"^-^ibj/bi;p) 



'rn-l' . 

n- 

n=0 



(4.18) 

where the prime on the product means as always that we omit the divergent term. Let us 
stress here that although in the analysis of this section we have used the jargon of theories 
of class S the final result (4.18) applies to any theory of the setup of figure 2. 



5. Properties of the operators &{r,s) 

We now proceed to investigate the algebraic properties of the difference operators ©(r,s) ob- 
tained in the previous section. These operators have many beautiful mathematical properties. 
Some of these properties are to be anticipated on purely physical grounds following the fact 
that &(r,s) compute resiudes for theories enjoying generalized S-duality. 
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We will find further evidence for our main claim, that ©(^.s) ^^iTfi?] is the index of a Tir 
in the presence of two surface defects along the two orthogonal circles of fixed by ji2 and 
j34. We will see that the two defects are naturally associated to the spin r/2 and spin s/2 
representations of SU{2). 



Factorization 

We claim 



(r,s) 



Pi 



®(r,0)6(0,s) 



pq\ 

t 



®(0,s)®(r,0) 



To show this we write 



®(0,s)6{r,0)^ 



2 ™*p2' 



n 



■jii-i 
n 

.71=0 



e{p^+'^j-^Hhi/hj]q) 



n 

. r?i=0 



0(^pnJ-n,^m+ra,-ra,^^.|l^..p^ 



nj—rii „m— m,- 



>j/bi;p) 



and then using 



^(x;g) = (-irx"(? 



(5.1) 



(5.2) 



(5.3) 



we establish (5.1). This factorization property has a natural physical interpretation. We have 
claimed that the operators S(r,s) introduce two surface defects, labeled by integers r and s, 
along the two orthogonal planes, and thus can be introduced (almost) without interfering 
with each other. The two defects feel each other's presence only through the proportionality 
factor of {^Y^ 1 which as explained in section 3.3 amounts to a redefinition of the IR r-charge. 



Commutativity of the operators 

The two operators ©(i.o) and 6(o,i) commute, 

[6(1,0), 6(0,1)] =0. (5.4) 

Mathematically this follows from simple theta- function identities. For example considering 
the term proportopnal to I(q2p2 a, . . . ) in the commutator we get 

Oiy';i)Oiia~';P) e{la~';p)6{±a~';q) _ 
e{a?-q) e{pa?]p) e{a?;p) e{qa^;q) ' ^ ^ 

etc. In more generality one can show that all the operators &(r,s) commute with each other. 
Physically this result is expected since these operators compute residues in theories which 
enjoy S-duality and we can employ this duality to change the order in which the residues are 
taken, see figure 4. 
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Figure 4: By generalized S-duality, the order in which one extracts the residues in fugacities a and 

b is expected to be immaterial, indeed the three different decompositions of the surface shown in the 
figure are topologically equivalent. This implies [&(r,s) ) *3(r',s')] ~ ^• 



Self-adjointness 

The operator ©(1,0) is self-adjoint under the natural propagator measure A • Xy, 

^ ^ A(a)Xy(a) /(a) [6(0,1) • 9{a)] = j A(a)Iv(a) [©(0,1) • /(a)] g{a) . (5.6) 

This can be shown simply by a change of the integration variable a ^ ag^5, if one assumes 
that the test functions /(a) and g{a) do not have poles in the strip \q\^ < \a\ < \q~^\- Under 
this assumption we write 



da 
a 

db 



Iv{a)A{a) f{a) 



-1/2^ 



(5.7) 



e{b-^q~^;p)- 



The integrand of the first integral satisfies 

Iyibq-^)A(bq-^) = ^(1 - qb^^l - b\~')PE 



qp-t 



b- 



PE 



1v{h)A{h) 



pt 



1-p 

e{b-^;p) 



PE 



(- 



p 



+ 



m _ + 

t ^ 



l-p l-q (l-p)(l-g) 



1-p 
){q-H'^ + qb-'^ - I) 



(5.8) 



and analogously for the second term thus establishing the self-adjointness of ©(0,1) • This 
property is to be expected from how we introduced the difference operator, namely as a way 
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Figure 5: Three different but equivalent ways to introduce a surface defect in a tube. First by 
coupling the two edges of the tube to a free hypermultiplet and computing a U{l)f residue of the 
index. Then by acting with the relevant difference operator on either of the edges and then gauging. 
The equivalence of these different procedures implies that the index is self-adjoint under the natural 
measure. 

to evaluate residues in the SU{2)f flavor fugacity of Tuv- Consider the setup of the first row 
in figure 5, where two theories T and T' are connected to a bifundamental hypermultiplet 
by gauging two SU (2) groups. The index can be written as double-contour integral for the 
two SU{2) gauge fugacities. To extract the residue in the SU{2)f fugacity, we can follow the 
procedure of the previous subsection and pinch one of the two contours for fixed values of the 
second integration variable. The result is the second or third column in figure 5, according 
to the order in which we perform the countour integrals. Since the setup is symmetric under 
exchange of T and T' , it follows that '3^r,s) must be self-adjont under the propagator measure. 



Acting w^ith difference operators on the trinion 

In the 2d TQFT interpretation of the index of theories of class S, the action of the difference 
operator ©(r,s) on a flavor fugacity corresponds to the fusion of the "degenerate puncture" 
associated to (r, s) surface defect with a flavor puncture. The choice of flavor puncture on 
which the operator acts is then expected to be immaterial. This follows directly from our 
construction of the operators and from the generalized S-duality enjoyed by theories of class S. 
The difference operators compute a residue of the index at a pole of one of the flavor fugacities. 
To obtain the operators we have singled out an additional flavor fugacity by decoupling a free 
hypermultiplet as in figure 3. Invariance of the index under S-duality means that our choice 
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of the additional flavor fugacity must be immaterial. The crucial claim is that the difference 
operators must give the same result when acting on any of the three legs of the elementary 
"trinions". By decomposing the Riemann surface into pairs of pants, it follows from this 
claim that the difference operators can be freely moved around the surface, indeed we have 
already shown that they are are self-adjoint with respect to the propagator measure. 

In the Ai case the index of the basic building block, the index of the trinion, is explicitly 
known, so we can establish this claim by direct computation. The Ai trinion index is the 
index of a free tri- fundamental hypermultiplet, 

I{a,b,c) = r (t^a^^b'^^c'^^;p,q) . (5.9) 



Dropping overall flavor-independent factors we obtain 



'e{\a-^-p)e{^\ab^^c^^-p) e{y;p)e{^la-^b^^c^^;p) 



(5.10) 

In the Macdonald limit p = [23] it is easy to show that the above expression is symmetric 
under permutations of flavor fugacities, 

©(0,1) (a)X(a, 6, c)~r(^y^a±i5±ic±i; 0,(7^ x (5.11) 
1 + - iXadjia) + Xadjib) + Xadjic)) [ J + (^) ] + («)x/ (&)x/ (c) (£j 



3/2 



We can also easily check the a ■v^ b ■v^ c symmetry Schur limit q = t [23, 22], 

9{a-^;p)e{ab^^c^\p) _ e{a'-p)9{a-^b^^c^^;p) _ j q,^,^^^^^,. . .5 

r{a^'b^'c^';p,q) = 0(„^±i^±i.p)^(„-i5±i,±i. ^) = yfpfpf^ify^)|pfpTpf7^ ' 
we obtain 

6(o,)(a)X(a,5,c)^^^^^^^^^^^^^_^^. (5.13) 

This expression is manifestly symmetric and also independent of p as expected [23]. For ar- 
bitrary (p, q, t) one has to verify that the combination of theta functions appearing in (5.11) is 



^The same expression can be also obtained thus from (5.11) by setting t = q. 
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symmetric under permutations of the three fugacities, an exercise we leave to our enthusiastic 
readers. We have checked this claim perturbatively in p. 

Even if the theory has no flavor punctures, we can still define an operation that modifies 
its index to include surface defects. We consider any pair-of-pants decomposition of the 
Riemann surface, act with the difference operator on any of the open punctures and glue the 
pieces back together. The result is well-defined, since the properties just discussed guarantee 
that it does not depend on the specific decomposition or choice of open puncture. 

In a 4:d theory with surface defects, it is generally difficult to separate the 4d and 2d 
degrees of freedom counted by the index. In the simple case of the theory of a free hypermul- 
tiplet endowed with surface defect, we can however identify an interesting operator intrinsic 
to the defect. Let us consider the "Coulomb limit" [23]^'' (see appendix A): t,p,q ^ with 
p'^/t and q'^/t fixed. The index of the trinion acted upon by (5.10) becomes 

©(o,i)(a)X(a,6,c)~ 1 + ^. (5.14) 

The non-trivial term comes entirely from the difference operator itself, and we identify it with 
a Coulomb branch operator X living on the defect, obeying the constraint = (as natural 
for SU{2) gauge group). Its 2d charges (see appendix A) are Lq = ^, Jq = —1, and / = 0. 
It is a naturally interpreted as the exactly marginal twisted chiral field that couples to the 
relative gauge coupling, which we mentioned at the end of section 2.2?^ 



Recursion relations 

Finally, the difference operators satisfy interesting recursion relations. Consider the computa- 
tion of S(o,s)®(o,i)^- Let Z[Tuv] be the index of the theory with one extra U{1) j puncture. To 
compute S(o,s)®(o,i)2^ extract the residue of S(o,i)X[7(/v'] at a = t^l'^q^l'^ . If we substitute 
the explicit expression for 6(o,i) we get the neat relation. 

This recursion relation can be formally recast as a fusion product between operators corre- 
sponding to representation of SU{2). We can write (5.15) as 

[s] X [1] = [s + 1] + [s - 1] , (5.16) 



^° A related limit for A/" = 4 SYM was discussed in [57] . 

■^^For higher rank theories there is an analogous story. For instance the Coulomb index of a bi-fundamental 
free hypermultiplet with a basic surface defect is X^flJo^ (^)*- Assuming that the Coulomb limit commutes 
with the residue computation this will be true for all indices. 



- 29 - 



where [s] is an operator corersponding to an s + 1-dimensional irrep of SU{2). This is 
reminiscent of the fusion relation between degenerate operator in Liouville theory, which in 
the AGT correspondence are associated to surface defects [58]. 

One can check (5.15) by using the exphcit expression for ©(o,s) (4.18). One then finds 
the same theta-function identity that we encountered in verifying the symmetry under per- 
mutation of flavor fugacities of the trinion when acted upon with a difference operator (5.10). 
This identity takes the following general form 

As suggested in section 3 it is natural to strip off a factor 

e{t;p) e{tq;p) e{tq'~';p) 

e{q~\p)e{q~^-p)"' e{q-^-p) 

to define an operator S(o,s)- Then the recursion relation takes a slightly simpler form, 



(5.18) 



©(o,.)6(o,i) - ©(o,.+i) + -J^Tp^e^uf^fio,s-i) ■ (5.19) 

Remarkably, the ratios of ^-functions can be interpreted as 2d indices X2d of free chiral fields 
with a specific R-charge assignment (see appendix A) 



e{^ix-p) 



X2d{x) . (5.20) 



e{x-p) 

In terms of these 2d indices we can write a final form of the recursion relation, 

X2d{t) 6(o,s)(3(o,i) = X2d{t q')&[Q,s+i) + X2d{t~^q^')&(Q,s-i) ■ (5-21) 



6. Bootstrapping the index for Ai theories 

We are now ready to implement the "bootstrap" of the 2d TQFT that computes the index 
of class S theories of type Ai. In other terms, we will determine the index from consistency 
conditions alone. Our strategy is a simplified (topological) version of the "Teschner trick" 
[24]. A degenerate puncture associated to a surface defect can be fused with any of the 
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flavor punctures, and assuming 5-duality the result must be independent of the chosen flavor 
puncture. 

To proceed, we make a plausible technical assumption: the difference operators ©(r,s) 
admit a complete set of eigenfunctions {ipx{a)}, normalizable under the propagator measure, 
with non-degenerate eigenvalues E^^,^y The label A runs over the irreducible SU{2) rep- 
resentations. Since the difference operators are self-adjoint in the propagator measure, an 
implication of this assumption is that any two different eigenfunctions are orthogonal in the 
propagator measure, and can be normalized to be orthonormal. This assumption is certainly 
true in the degenerations limits p = or g = 0, where the difference operators become the 
familiar Macdonald operators^^ and the eigenfunctions are proportional to Macdonald poly- 
nomials. For arbitrary (p, q, t), to the best of our knowledge the eigenfunctions are not known 
in closed form: they are expected to be an "elliptic" deformation of the Macdonald functions. 
In appendix B we describe an approximation scheme to determine the eigenfunctions in an 
expansion for small q and p. 

Expanding the index of the trinion, 

Xo,3 = Yl ria)^^ibW{c) , (6.1) 

and acting with the difference operator one any of the flavor punctures, we have 

As by assumption the eigenvalues are non-degenerate, this implies that the index is diagonal 
in the basis {tp°^}, 

:io,3 = 2] r (a)r (^)r (c) . (6.3) 

a 

It remains to fix the structure constants Ca- Using orthonormality of the eigenfuctions under 
the propagator measure, the index of the theory associated to the four-punctured sphere is 
immediately given by 

= Y.^l r{a)r{b)r{c)r{d) . (6.4) 

a 

^^Of course, for the general Ai theory, which has a Lagrangian description, the index is already explicitly 
known in terms of a matrix integral. Since everything is explicit, S-duality of the index for the Ai quivers can 
be rigorously established [21] using an identity for a certain integral of elliptic Gamma functions [59]. However, 
the bootstrap approach generalizes easily to higher-rank theories, whose index is a priori unknown, and even 
for the Ai theories it has the virtues of directly giving the index in the very useful "diagonal" form. 

■^^More precisely, they are related to the Macdonald operators by a similarity transformation. 
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According to our prescription, extracting the residue at a = t2 closes a puncture, so 

1 „ 



=-0,3 



We conclude that 



2 Iv Res„^ji/2 - xo,4 
a 



1 



(6.5) 



2Iv (Res,=,v. ir(a)) ' 

The index of general Ai theory associated to a surface with k flavor punctures is given by 

k 



(6.6) 



(6.7) 



i=l 



To introduce an (r, s) surface defect, we are instructed to simply multiply by E'^j. inside the 
summation sign. 

We can also give expressions for the eigenvalues using the residue interpretation of the 
difference operators. For example 



171 a 



Res 11, ^V'"(a) 



Res 1 

a=f2 



(6.8) 



In particular in the Macdonald limit, p = 0, we have 
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(to±2;g) 



.i=0 
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l-t'_q_ 



P"(t2g2^;g,i) 



(6.9) 



Here P°'(z] q, t) are Macdonald polynomials (B.5) and we have defined as always 

00 



(6.10) 



^=0 



It is easy to verify that then the energies satisfy the recursion relation 



^(0,s) ^(0,1) 



1-q 



1 - qH 



-E, 



(0,s+l) 



+ 



l-t 



1 - q-H- 



TP Q 



(6.11) 



which is a direct consequence of recursion relation satisfied by the operators (5.15). 

In the Macdonald limit the energies £"(^,0)" can be related to modular S-matrix of the 
refined Chern-Simons theory [60]. Up to a independent normalization factor these energies 
are given by 



P''it-2q-2';q,t) S^' So 



(6.12) 
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where Sa"^ is the modular S-matrix of the refined Chern-Simons theory as defined in [60]. In 
the Schur hmit {q = t) the refined Chern-Simons theory further reduces to the 2d q-deformed 
Yang-Mills (in the zero-area limit), as defined in [61]. The relation of 2d qYM to the index 
was discussed in [22, 23]. In the 2d YM language the flavor punctures correspond to fixing the 
holonomies of the gauge fields around the punctures. Thus, at least in the Macdonald limit, 
introducing a surface defect corresponds to adding a puncture with fixed holonomy for dual 
variables of the 2d theory, i.e. the canonical momenta dual to the gauge fields [62, 63, 61]. 



7. 3d reduction 

Starting from the Ad index, which is the supersymmetric partition function on 5^ x 5^, we can 
dimensionally reduce on the to obtain the partition function of a 3d theory on a squashed 
. The reduction of the 4d index to 3d partition function is achieved by the following scaling 
of the fugacities, 

q = e'^'^\ p = e'^'^\ t = e'^ , = e^'^'"^ (7.1) 

and taking /3 — )• [64, 65, 66], where ai are 4d flavor fugacities and mi 3d real masses. The 
4d index goes to the "ellipsoid partition function", with squashing parameter 

5=^^. (7.2) 

UJ2 

The 3d theory which arise in the infrared from the circle reduction a id theory of class S has 
a nice mirror description [67] in terms takes of a star-shaped quiver: a central node connected 
to a set of linear quivers, where each linear quiver is associated to a puncture of the Riemann 
surface. The wavefunction 'ip°'{a) that we associate to a puncture in the TQFT description 
of 4d index reduces in the 3d limit to the partition function of the corresponding linear 
quiver tail [68]^^. 

The surface operators of the Ad theory are expected to become line defects of the 3d 
theory. Following the chain of dualities, one finds that the canonical surface defects for a 
4d theory of class S, which arise from codimension four defects of the 6d (2,0) theory, go to 
Wilson loops for the central node of the quiver in the mirror 3d description. We would like 
to verify this fact. 

For Ai theories the mirror description involves a 3d gauge theory with a U{1) gauge 
group for each puncture, and a single SU{2) gauge group. The matter content includes 
a doublet of hypermultiplets for each puncture, with charge one under the corresponding 



See also [691 for related work. 
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C/(l) gauge group, and g adjoint hypermultiplets for the SU{2) gauge group, where g is the 
genus of the Riemann surface. Each "quiver tail" consists of a U{1) gauge theory coupled 
to a doublet of hypermultiplets. This is the so-called T[SU{2)] theory, which is intimately 
related to the S-duality operation in four-dimensional = 4 SYM [70]. In particular, its 

partition function is related by a generalization of AGT [71] to the integral kernel which 
implements the S operation on conformal blocks on a one-punctured torus. This fact has 
a simple consequence: the partition function of T[SU{2)] satisfies a functional equation, 
and intertwines the action of two operators, which are related to the action of 't Hooft and 
Wilson loops of = 4 SYM. 

The Wilson loop operator multiplies the T[SU (2)] partition function by a character of 
the SU(2) flavor symmetry which rotates the doublet of hypermultiplets. In particular, the 

partition function of the full star-shaped quiver in the presence of a Wilson loop operator 
for the central SU{2) node can be obtained by acting with the Wilson loop operator on any of 
the T[SU{2)] linear quiver tails partition functions. This can be then rewritten as the action 
of a 't Hooft line operator on the partition function of that quiver tail. We will now show that 
3d limit of the operator which inserts a surface defect in the 4(i theory exactly reproduces 
the 't Hooft line operator. This will verify that the insertion of the surface defect goes to the 
insertion of a Wilson line defect for the middle node of the mirror star-shaper quiver. 
Let us consider the 3d limit (7.1) for the operator ©(o,i)- The prefactors 



e{a^-p) sin27r^ 



(7.3) 



while shift operators reduce to 



f{q^a) ^ /(m+^a;i). (7.4) 



All in all. 



sin^ ^-^l"^'" sin^ ^-^^+^'" 



^ sin sin27r^ A„_,^„i^^ . (7.5) 

(In complete analogy, the limit of the operator ^[ifi) is obtained by switching uji with uj2.) 
Under an appropriate identification of the parameters this coincides with the 't Hooft line 
defect, as desired. Line operators in the framework of AGT were studied in [72] for the Ai 
case and in [73] for higher rank theories. One can compare the operator (7.5) to equation 
(5.25) in [72], which displays operator inserting a 't Hooft operator in the J\f = 2* theory. 
The ratio of the two periods is the b parameter of the underlying Liouville CFT, and the 
modular properties we will discuss below are related to the change 6 — )• 

Finally let us briefly discuss some mathematical properties of the 3d difference operators. 
The 3d reduction of g) is nothing but the basic Macdonald difference operator, with 
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effective parameters 



t = e "2, q = e '^2, a = e '^2. ( ' -oj 

Similarly, for the reduction of 6(0,1) we get the Macdonald operator with parameters 

f = e "1 , q = e '^1 , a = e "1 . (7.7) 

Note that q and q' are related by a "modular" transformation, 

g = e2^^ ^ g' = e-^. (7.8) 

Physically this follows from the fact that the (1,0) and (0,1) surface defects have support on 
the two maximal (linked) circles in fixed respectively by ji2 and j^^. Interchanging the 
two types of defects is the same as interchanging the circles. Let us also write down the 3d 
reduction of the elliptic measure, 

^ 2 r(a±2;g,p) ^ 2 17,+ij+i-l±2ma+ii+i±2m' ^'^ 
where we have defined Qij = iui + juj2 ■ 



8. Aat-i theories 

In this section we outline the generalization of our "topological bootstrap" procedure to 
the higher rank theories. We first discuss, but not review in detail, the six-dimensional 
construction of = 2 four-dimensional gauge theories and surface defects in them, and the 
correspondence with two-dimensional Toda CFT correlation functions. We refer the reader 
to [16, 74, 75] for a complete discussion of the relevant facts. 

The (2, 0) Qd SCFTs are labeled by a simply-laced Lie algebra g. An important class of 
codimension-two (to be contrasted with codimension-four) superconformal defects is labeled 
by an embedding p of su{2) in g, see e.g. [74, 76]. The properties of these defects can be 
motivated, for example, upon compactification on a circle shared by all the defects. The 
bulk theory becomes bd SYM with gauge algebra g. The codimension two defect induces a 
singularity of the bd SYM fields which breaks the gauge symmetry to g{p) (possibly times 
some Abelian factors) at the defect. The codimension-four defect becomes a Wilson loop 
in the representation of 5 or g{p). In the construction of four-dimensional M = 2 gauge 
theories, one considers the twisted compactification of the Qd SCFT on a Riemann surface C, 
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with codimension two defects at points in C. Codimension four defects at points in C produce 
surface defects in the 4d theory. 

Nekrasov's instanton partition function for the 4(i theory is conjecturally the same as a 
conformal block for the Wn current algebra [20, 77]. The vertex operators for the Wn theory 
come in families associated to various patterns of degenerate vectors. The families are labeled 
by the embedding p of su{2) in g, and an appropriate choice of Toda momentum, a vector 
valued in the dual of the Lie algebra g. Non-degenerate vertex operators (corresponding to 
/? = 0) are labeled by an imaginary Toda momentum. Semi-degenerate vertex operators are 
labeled by the sum of an imaginary momentum in f{p), and a real momentum of the form 
bX + b~^X' , where A, A' are weights for g{p)- For a fully degenerate vertex operator, f{p) = 
and g{p) = g. 

The vertex operators map to a six-dimensional configuration with a codimension two 
defect labeled by p and codimension four defects labeled by A and A', sitting on the 12 and 34 
planes in flat space. In the Toda context one can change the type of a defect, by analytically 
continuing the continuous imaginary part of the momentum to some discrete choices of real 
momentum. 

Much of our analysis concerning the analytic structure of the index can be repeated for 
the analytic structure of the Nekrasov's partition function, i. e. for the properly normalized 
conformal blocks. The conformal blocks have poles at discrete values of the Toda momenta, 
which correspond to special values of the mass parameters of the four-dimensional theory. 
At these values, the divergence is caused by a flat directions in the Higgs branch. Thus 
the process of making a Wjsr vertex operator more degenerate can be matched to a Higgsing 
process, in the presence of vortices which become surface defects in the IR. 

Let us start from the basic example involving a bifundamental hypermultiplet. Our con- 
struction starts from a theory Tir, which for simplicity we can take to have a six-dimensional 
description in terms of the ^at-i theory on a Riemann surface C with at least one "full" 
puncture, which carries the SU{N) flavor symmetry. Then Tuv is a theory associated 
to C with one extra simple puncture inserted near the full puncture. The position of the 
simple puncture controls the gauge coupling of the auxiliary SU (N) gauge group. The simple 
puncture has a U{1) flavor symmetry, which we gauge in the Higgsing procedure. 

In the context of conformal blocks, the Toda momentum of a simple puncture is an 
imaginary multiple of a single simple weight wi. It can be analytically continued to the real 
momentum of a fully degenerate puncture, but only one with X = rwi and A' = swi. These 
surface defects are naturally associated to r-th and s-th symmetric powers of the fundamental 

^^This assumption does not remove any generality from om analysis: no matter what Tir is, one can build 
an auxiliary six-dimensional theory in terms of a sphere with two full punctures, and Tir added by hand at 
one puncture, by gauging the diagonal SU{N) flavor symmetry. 
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representation. In the context of the superconformal index, they are the defects introduces 
by the difference operators &(^r,s)- 

Can we discuss surface defects labeled by other representations of SU (N)? In the context 
of 5d SYM, Wilson loops in representations Ri and R2 can be brought together to a Wilson 
loop in the representation Ri <^ R2, which can then be decomposed to a sum of Wilson loops 
in irreducible representations 

Ri0R2 = ^Ni2Ri. (8.1) 

i 

If we try to uplift such a formula to an "OPE" of codimension-four defects in 6(i, 

SR,^SR, = Y,rl2 0SR^, (8.2) 

i 

we need to promote the coefficients to 2d theories 7^2 with Here "OPE" is intended as 
a collision of the punctures in C, not in space-time: the configuration of two codimension-four 
defects represents a single, two-parameter surface defect in four-dimension, which simplifies 
to a sum of defects if the parameters are tuned appropriately. 

At the level of the index, the location of the defects on C is immaterial, and the relation 
will become 

&{RuO) ® ©(^2,0) = Yl ^['^s] ® e(iJ„o) , (8.3) 

i 

where x['7i2] is the 2d index of the "OPE coefficient theory" . We have discussed such relations 
in the Ai case (see section 5). 



Closing minimal punctures 

In this section we discuss the generalization of the residue computation we performed for Ai 
quivers to higher rank cases. The generalization is straightforward and we will only outline 
the essential steps. As before, we consider a theory of class S Tir{= Tc) with at least one 
maximal puncture and couple to it a free bi-fundamental hypermultiplet. The index of the 
free bi-fundamental hypermultiplet is equal to 

Ihyp{h,c;a)= [[ [[ — ^— — ^. (8.4) 

Here hi and Cj (11^=1 — YYj=i'^j — 1) SU{N) fugacities and a is a U{1) fugacity. To 
couple this hypermultiplet to a general Aj^j^i quiver corresponding to Riemann surface C we 
gauge a diagonal symmetry SU {N) by integrating over fugacity c, 

Iuv= fll ^;^A{c)Iv{c)Ihyp{h,c-a)Ic{c-\...). (8.5) 
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Figure 6: To compute the residue of the index at a pole of the U{1) fugacity we consider the 
degeneration hmit of the theory by decouphng a free-hypermuhiplet. 

This index has poles at the following values of the U{1) flavor fugacity 

a = tl//^g^/^, (8.6) 

with non-negative integers r and s. Actually, there are poles also at a — )• exp[^^] a as we 
discussed before (see section 3.1). A way to see this is again by studying the pinchings of c 
contour integrals. Let us start by considering the poles of the integrand of (8.5) at 

c, = i^g™'p"»-^, i = l,...,N. (8.7) 

Here a{-) is an element of S{N). Notice that these combined poles have a straightfor- 
ward physical interpretation, as contributions from a region in field space where complex 
bi-fundamental fields Q^^*^ receive a vev, which Higgses the U{N) gauge symmetry com- 
pletely. The SU{N) gauge invariant operator is a determinant of the bifundamental fields. 
Thus the poles at a^t^/'^p^q^ = 1. The residue (multiplied by N for future convenience) is 
Xy^ Xfjj^ih, . . .). If we gauge the U{1) flavor symmetry with fugacity a, the contribution to 
the index of the simple pole at a = t2 is N'^^Zj-j^^ih, . . .). The other poles at a = exp[^^] ta 
give the answer A^~^l7-^^(exp[— b, . . . ), and they combine to give the Zj^j center gauging 
of TiR, as expected (see section 3.1). 

Further specifying to a = t^pNqN we get that the simple poles above pinch the inte- 
gration contour and become double poles on the unit circle, because we have only — 1 
independent SU (N) flavor fugacities, if 

N N 

s = '^mi, r = ^ni. (8.8) 

2=1 i=l 
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This implies that the index has simple poles at (8.6). The residue computation is given by 
the action of the following operator on the index 



'^[TlR,G[r,s)\ — ®{r,s) '^[Tir\ — 



Ei=i"»=^Ei=i"^i= 



TV 

n 



'nii-l 

n 

. m=0 



e{p-^jq^-^jhj/bi;p) 



Tli — l 

n 

.n=0 



The basic properties of these operators are the same as for the Ai case. These operators are 
all commuting and self-adjoint under the natural measure. 



The index of a general theory of class S 

We are finally ready to bootstrap the index of class 5 theories of type ^tv^i. The differ- 
ence operators 6(r,s) have common eigenfunctions {V'"(a)}, labelled by irreducible SU{N) 
representations. The eigenfunctions can be chosen to be orthonormal under the propagator 
measure and we assume again that their eigenvalues are non-degenerate. By the usual logic 
of introducing a degenerate puncture and colliding it with flavor punctures, we deduce that 
the index of the trinion with three maximal punctures is diagonal in the {ip°'{a)} basis, 

Itrinih, C, d) = E Ca V'"(b)V'°(c) V°(d) . (8.10) 

a 

The determination of the structure constants is somewhat more involved than in the Ai case. 
Let us outline the procedure for the case of A2 ■ 

The index of the free hypermultiplet of A2 can be expanded as 

X,,,p(b,c;a) = J^V"(b)r(c)(/>f2,i)(a). (8.11) 
a 

Here b and c are SU{3) fugacities and a is U{1) fugacity. The label (2, 1) denotes the 
auxiliary Young diagram corresponding to the U{1) puncture. This equation defines the 
function (p^^i-^ia) since the left-hand side is explicitly known (8.4). (That Ihyp(h,c;a) must 
take this diagonal form follows again by consistency of the theory with an extra degenerate 
puncture.) Consider now the theory associated to a sphere with two maximal and two minimal 
punctures. The index can be calculated in two different duality frames. In the degeneration 
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limit of the surface where a maximal puncture collides with a minimal puncture, the theory 
is obtained by gauging the diagonal SU{3) of two bifundamental hypermultiplets, 

2 

^=E'^(2,i)(^/^)<^f2,i)((c&)"') n^"(^^) ■ (8.12) 

a 1=1 

Here h and c denote two U{1) fugacities. On the other hand, in the degeneration limit where 
two maximal punctures collide, the trinion with three maximal punctures is coupled to a free 
hypermultiplet by gauging an SU{2) subgroup (this is the duality discovered in [78]). In this 
duality frame the index is given by 

^ I' dz 

^ = E^" n^"(^^) f ^A(z)Xy(z)X;,,p(z,63)V;°(zc,Z^^C,C-2)_ (8_^3) 

a e=i 

Here ihyp{z, b^) is a free hypermultiplet in the fundamental of SU{2), where the SU{2) 
fugacity is parametrized by z\ the two half-hypermultiplets have charges ±3 under the U{1) 
symmetry parametrized by h. The U{1) charges of the different components in the two duality 
frames are linearly dependent [78], which results e.g. in the non-trivial power of h in the free 
hypermultiplet Xhyp. Comparing the indices in the two duality frames we obtain 



'^f2,l)(V^)'^f2.1)((^^)~'^ 

§ ^^A{z)Iv{z) ihypiz, b3)i;»{z c, z-i c, c-2) 



^= . (8-14) 



The right-hand side here is written as a function of b and c but is actually independent of 
these as an implication of S-duality; this can be explicitly checked [23]. 

Finally, by gluing elementary building blocks, we can write the index of a genus g theory 
with k maximal and k' minimal punctures, 

k k' 

J2iCa)''^'^'^'' n^"(^^) n ^i2,l)ibm). (8.15) 
a £=l m=l 

This algorithm can be easily generalized to higher rank. The trinion with three maximal 
punctures is always diagonal in the {^/'"(a)} basis, see (8.10). The index of the free hypermul- 
tiplet has a diagonal expression analogous to (8.11), in terms of two tp'^ wavefunctions and 
one wavefunction (j)^^ ^ associated to the minimal puncture, which can in fact be fixed by 
the known expression of the index of the free hyper. Finally the structure constants Ca are 
determined by comparing two degeneration limits of the N + 1-punctured sphere with two 
maximal and — 1 minimal punctures: in one duality frame we have a linear superconfor- 
mal quiver with the two maximal punctures at the ends, and in the other one we have the 
SU{N) trinion (the Tjv theory) coupled to the superconformal tail with a maximal puncture 
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on one end and minimal puncture on the other end. The two duahty frames are related by a 
generalized Argyres-Seiberg duality [15]. The former duality frame is completely Lagrangian 
and the index can be computed explicitly. The index in the latter frame is written in terms 
of the index of the trinion. Equating the indices in the two duality frames one determines the 
trinion structure constants C^. These date are sufficient to fix the index of any theory con- 
taining only maximal and minimal punctures. More general punctures can be incorporated 
by a Higgsing procedure that involves more general superconformal tails. In appendix C we 
give an example of this procedure. The final prescription is spelled out in [43] and coincides 
with the prescription of [23].^^ 



9. Discussion 

In summary, the superconformal index, being a protected quantity, must be consistent with a 
large class of deformations of the 4c? field theory. Moving along an exactly marginal direction 
the index does not change, so it must be the same when evaluated in different S-duality 
frames. This has the non-trivial implication that the index of class S theories is computed 
by the correlator of a 2d TQFT [21]. In this paper we have considered a more interesting 
class of deformations: we have studied how the index is affected by expectation values of 
supersymmetric operators. We have argued that a pole of the index in fiavor fugacity is 
associated to a bosonic fiat direction, parametrized by the vacuum expectation value of a 
protected operator, and that the residue captures the IR physics reached at the endpoint of 
the fiow triggered by the vev. We have formulated a precise prescription to evaluate the index 
of the IR theory. 

We have focussed on a special class of poles, associated to RG flows that introduce 
two-dimensional defects in the IR SCFT. In particular, starting from the index of a theory 
with fiavor symmetry Q x U{1) our residue calculus determines the index of the theory with 
smaller fiavor symmetry, Q, but endowed with BPS surface defects. In the language of the 
2d TQFT, the surface defects are associated to special "degenerate" punctures. The fusion 
of the degenerate punctures with the "ordinary" punctures associated to fiavor symmetries 
amounts to acting on the fiavor fugacities with certain elliptic difference operators. Under a 
minimal set of assumptions, namely consistency of the 2d TQFT structure and knowledge of 
the index of free theories, we were then able to determine the index of a general theory of 
class S. 

■^^The prescription of [23] gave divergent results for a small subclass of theories of class S, however following 
the residue logic described in this paper it was explained in [43] how to resolve those singularities. 
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We can perhaps extract two general lessons for quantum field theory. First, in a given 
theory, even if ultimately we are most interested in the familiar correlators of local operators, it 
is fruitful to consider the larger set of observables that includes defects of various codimension. 
In our problem the introduction of surface defects was the key step. Second, to obtain results 
in a given theory, it it fruitful to enlarge the view to all the theories of the same class. Indeed 
our "bootstrap" is a conventional bootstrap for the 2d TQFT, but for the Ad SCFT it is really 
a bootstrap in theory space. 

Some open problems 

In closing, let us briefly mention several further research directions, questions, and specula- 
tions arising from our work. 

• First and foremost, it would be very valuable to have a bona fide computation of the 
superconformal index for Lagrangian class S theories endowed with surface defects. We 
have in mind a first-principles calculation directly from path integral definition, using 
the methods of supersymmetric localization. 

• The difference operators that introduce surface defects are sums of terms shifting the 
values of flavor fugacities in different ways, with coefficients looking as indices of 2d free 
fields, i.e. combinations of theta-functions (see e.g. (8.9)). 

It is natural to wonder this sum has a more direct physical interpretation. What happens 
physically when the defect puncture collides with a flavor puncture? There is some 
evidence [75] that the surface defect may degenerate to the direct sum of simpler, 
decoupled defects. Schematically, we could write 



where at are the simpler surface defects, and the "coefficients" Tj are further decoupled 
2d degrees of freedom. We can imagine that the simple defects cjj give the individual shift 
operators in the index fugacities, while the 2d index of the Tr^i produces the coefficients 
in the expansion of the difference operator. It would be interesting to make sense of 
this conjectural structure. 

Similarly, we have seen that the difference operators introducing the defects satisfy 
recursion relations. The coefficients of these recursion relations are again combinations 
of theta-functions, and could perhaps be interpreted as indices of 2d theories. 

• The index of class S theories takes a diagonal form in the eigenfunctions ip°'{a) of 
the *3(r,s) difference operators. The eigenfunctions do not depend on the specific four- 
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dimensional theory. It is tempting to interpret these wavefunctions as universal modules 
of some intricate algebra, built out of the flavour currents and their descendants. 

• Finally, it will be important to generalize our results N = 1 theories. The simplest 
route is to take the N = 1 ("Sicilian") limit of our construction, which corresponds to 
giving a mass to the chiral field in the N = 2 vector multiplet [79]. At the level of the 
index this amounts to setting [80] pq = t^. It would be interesting to check if the index 
of surface defects remains sensible after this mass deformation to A/" = 1. 
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A. Embedding 2d in Ad 

Let us discuss the embedding of 2d (2, 2) superalgebra in the M = 2 Ad superalgebra. We 
choose to parametrize 5^ x by = 1 for S'^ and by an xq for the . We start 

from the J\f = 2 superconformal algebra S'C/(2,2|2), whose bosonic subgroup is SU{2,2) x 
SU{2)r X U{l)r. A maximal subgroup of 5f/(2, 2) is SU{2)j, x SU{2)j^ x U{1)e. The 
generators ji2 = j'l + j2 is the rotation of the (xi, X2) plane and ^34 = j2 — ji is the rotation 
on the (2:3, 3:4) plane. 
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Table 1: For each supercharge Q, we hst its quantum numbers, the associated 6 = 2 {Q, Q^}, and 
the other 5s commuting with it. Here / = 1, 2 are SU (2)r indices and a = ±, d = ± Lorentz indices. 
E is the conformal dimension, (ji,j2) the Cartan generators of the SU{2)i SU{2)2 isometry group, 
and {R,r), the Cartan generators of the SU{2)ji (g) U{l)r R-symmetry group. 

The two dimensional defects wrap the and one of the equators: either xi = X2 = 
or X3 = X4 = 0. The surface defects preserve an SU{1, 1|1) x SU{1, 1|1) x U{l)j subalgebra, 
where the U{l)f commutes with aU the preserved supercharges and can be considered as a 
flavor symmetry from the 2d perspective. The bosonic subalgebra is given by two copies, left 
and right, of SU{1, 1) x U{l)j and by the flavor [/(!)/. The SU{1, 1) is a chiral half of the 
global conformal group in two dimensions with Cartan generator Lq (^o); while the U{l)j 
is generated by Jq (Jq). We denote the generator of the flavor ^7(1)/ by /. The fermionic 
generators of the two SU{1, 1|1) NSNS algebras are denoted by G^, Gf with r = 

There are many equivalent ways to embed this into the 4(i algebra but since we are 
computing the index with respect to Q^^^ there are two choices which correspond to two 
possible choices of /, 

f = ji2 + R, f' = m + R. (A.i) 

The two choices are of course isomorphic and correspond to putting the defects on the two 
different equators. For concreteness we will focus on the former choice. The left and the right 
super-charges are then mapped as follows 

(G+i , G: = (Qi-, , (G+ 1 , G: = (Qi^, Q2+) . (A.2) 

22 22 

The Gt and Gf are the superconformal counterparts. The map between charges 
2 2 

Jo = i? + r + /, J^ = R-r + f, (A.3) 

Lq + Lq = E , Lq — Lq= j34 , 
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generates the canonical (2, 2) superalgebra 

[Lo, Gf] = -rGf, [Lo, Gf] = -rGf, (A.4) 

[Jo,Gf] = ±Gf, [Jo,Gf] = ±Gf, 

{Gt , J = Lo ± J Jo , {Gt , J = Lo ± J Jo . 

2 2 ^ 2 2 Z 

The admixture of the flavor symmetry / inside Jo and Jo is fixed by the last anticommutators. 

Let us now consider our Ad index, 

I = Tr(-l) V"'' g^'*"'' f"^^ ■ (A.5) 
We notice that we can identify the charges as, 

j34 - r = Lo - ^Jo - (Lo - ^Jo), ji2-r = 2f - Jo r + R = Jq - f . (A. 6) 

Moreover only states with Lq = ^ Jo contribute to the index due to (3.3). Thus, the index 
can be written as follows 

X = Tr(-l)V"+^^"^-^P^ f — V ^ • (A.7) 

\pqJ 

This is an NSNS index of the (2,2) theory. Choosing the symmetry /' in (A.l) to be the 
flavor symmetry the index becomes 

I = Tr(-l) V°+^^''^-^' / ( -) , (A.8) 

\pqJ 

i.e. we just switch between p and q, and / with /'. 

We can give a 2d interpretation of the prefactor that we obtained in the residue compu- 
tation, 

where u is non-negative integer. The numerator can be understood as the partition function 
generated by two fermions with the following charges 

^ : Lo = -^u, f = -l-u, Jo = -2-n, (A.IO) 

tp" : Lo = l + ^^x, f = l + u, Jo = 2 + u, (A.ll) 

and their derivatives. The charges of the derivative which contribute to index are Lo = 1, 
Jo = 0, / = 0, and Jo = Lo = 0. This derivative contributes a factor of q to (A. 7) and 
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generates the infinite product in the theta-functions of (A. 9). Similarly, the denominator is 
generated by bosons with the following charges 

p-i-" : Lq = -1-Iu, f = -l-u, Jo = -l-u, (A.12) 

: Lo = ^ + ^u, f = l + u, Jo = l + u. (A.13) 

These are unusual charge assignments. The standard charges of a free chiral field (free complex 
scalar and a complex fermions) are obtained in the formal limit u = —1: here the flavor 
charges are zero and the expression (A. 9) diverges and hence the limit is formal. We take this 
as further indication that these partition functions describe decoupled fields with accidental 
symmetries, and that they should be stripped off from the definition of the surface defects. 
If we redefine the charges as 

Lo = Lo-^/, Jo = Jo-/, (A.14) 

then the new charges have natural values for the fields contributing to (A. 9). 

Note also that the combinations of theta-functions appearing in the basic difference oper- 
ator (4.16) are 6^^; q)/6{b^'^; q) which is again a free chiral field, with standard charges, 
also charged under fiavor symmetry coupled to b. 

Finally, we consider different supersymmetric limits of the index. For concreteness we 
consider the index (A. 7). 

"Higgs" index 

We take the limit g — )• while keeping the other fugacities finite. This is the 2d version of 
the Ad Macdonald index [23]. The index becomes 

Ih = Trn (-1) V^""^" t-^°~^ ■ (A.15) 
Here Tth is summing over states satsfying 

Lo = ^Jo, U = \jo. (A.16) 

Thus this index is half-BPS and gets contributions from states annihilated by two super- 
charges, one left and one right, of same R-symmetry charge. The index of the free 2d chiral 
field (A. 9) here becomes 

For instance the index of the Ai trinion with a surface defect in this limit is given in (5.11) 
and gets contributions from states charged under flavor symmetry and thus captures Higgs 
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branch states. Note that the 2d (2, 2) algebra does not constrain the value of charge / and 
thus there is no analogue here of the Hall-Littlewood index of [23] . 



"Coulomb" index 

We take the limit {q,p,t} — )• while keeping p/q and pq/t finite. This is the 2d version of 
the 4(i Coulomb index [23]. The index becomes 

Here Trc is summing over states satsfying 

Lo = -^Jo, U = \jo- (A.19) 

Thus this index is half-BPS and gets contributions from states annihilated by two super- 
charges, one left and one right, of opposite R-symmetry charge. The index of the free 2d 
chiral field (A. 9) here becomes 



Oitp"', q) (pq 



l+u 



1 _ M / 2 



(A.20) 



The index of the A\ trinion with basic surface defect in this limit (5.14) gets contributions 
from flavor singlets not counted by the Higgs index and thus naturally captures Coulomb 
branch physics. 



"Schur" index 

We take the limit p = t while keeping all the fugacities finite. This is the 2d version of the 
4d Schur index [23]. The index becomes 

= TV (-l)^g^o-5-^>-^ (A.21) 

Since / commutes with all the charges in the (2, 2) algebra this index is actually independent 
of q. In other words this index is simply Tr( — 1)^ refined only by flavor symmetries from 
the 2d perspective. Given such an index care has to be taken in expanding it in fugacities 
since there is no a-priori good expansion parameter. It gets contributions only from states 
satisfying 

Lo = Ijo, U = ]^Jq. (A.22) 
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The index of the free 2d chiral field here (A. 9) becomes 



, (A.23) 



and in particular all dependence on q factors out automatically as it should. Note that here 
the vacuum canceled out against a fermionic state with zero /-charge: this is an artifact of 
the issue mentioned above. 



B. Wavefunctions of the Ai elliptic RS model 

The difference operators that we defined through our residue prescription are closely related 
to well-known difference operators. Consider the similarity transformation 



Oit;p) 

where H is given by 



&^o,i^ = -^pl^^JCib)HIC-\b), K:{b)=r{tb^^;p,q) , (B.l) 



eitb-^;p) 

This operator is known as the elliptic Ruijsenaars- Schneider (RS) difference operator, which 
is a relativistic generalization of the Hamiltonian of the elliptic Calogero-Moser-Sutherland 
model [31, 32, 33]. In fact we have a second set of operators obtained by interchanging p 
with q. All these operators are commuting as we have discussed in the bulk of the paper. 
The spectrum of the RS operator is well-known in certain degeneration limits. Taking p = 
the operator reduces to the Macdonald difference operator, whose eigenfunctions are the 
Macdonald polynomials. Conjugating with fC the propagator measure (the Haar measure 
times the index of the N = 2 vector multiplet), we find 

.J ^1 \ ^ i r(t5^^;p, q) 

""^^ 2T{tb^^-p,q)T{b±^-p,q) ^ 2 T {b^^- p, q) ' ^""'^^ 

The operator H is self-adjoint under this measure. Setting p = this becomes 



^^A relation between Calogero-Moser-type models and elliptic Gamma functions, which are the building 
blocks of the index of the chiral superfield, was discussed in [81]. 
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which is the measure under which Macdonald polynomials are orthogonal. The Macdonald 
polynomials are explicitly given by 



i=X 



1 - gJ+i 
tqj 



(B.5) 



A i — 1 ^ -j A— X— 1.- A 

l — tQ-l -i-r I — tq-' 



i=0 j=0 



n 13^^^ 



2j-A 



It is important that these eigenfunctions are non-degenerate. 

For general {p,q,t), although there is some discussion of wavefunctions in the litera- 
ture [82, 83, 84] (see also [85, 86, 87]), to the best of our knowledge no explicit closed form is 
known. However, as a proof of concept that such eigenfunctions exist and are non-degenerate 
we can construct approximate wavefunction as follows. We solve for the eigenfunctions using 
the following ansatz 

a'+^g{a)-a~'"^g{a-^) 



a — a 



(B.6) 



Plugging this ansatz into the difference equation (B.2) we obtain. 



a^g{q^^a)0{ta?; p) a ^g^q^a ^)9{ta ^; p) 



+ 



(1 -g-ia-2)0(a2; p) (1 -g-ia2)0(a-2; p) 



+ 



(B.7) 



-q 2* 



a ^g{q 2 a ^)9{ta?]p) a^g{q 2a)6{ta 2; p) 



(1 - qa^)e{a^; p) 



+ 



(l-ga-2)^(a-2;p) 



Sr. 



a^~^^g{a) — a ^ ^9(0- ^] 



a — a 



As a first approximation we neglect the term in the first line, assuming as always that \q\ < 1, 
and obtain the simple equation 

go{q-^a)e{ta~'^; p) _ ^i^^f^Q^ go{a) 



(B. 



(1 -ga-2)0(a-2; p) ^ ^1 l-a-2- 

The subscript on goi') signifies that we are making an approximation here. This equation 
is solved using the following property of the elliptic Gamma function 



^iqz;p,q) = 9{z;p)T{z;p,q) . 



The solution is 



9o{a) = (1 - a ) 



2n r(a '^■,p,q) 



2*^ 



(B.9) 



(B.IO) 



r{ta '^■,p,q) ' 

This solution is a joint eigenfunction of the two hamiltonians ^^^d 6(0,1)5 ^-e- it is 

symmetric in p and q. The eigenfunction are given thus by 



, . ^ r(Q ^■,P,q) , r(a2;p,g) 



r{ta ^■,p,q) 



T{ta^;p,q) ' 



(B.ll) 
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Note that this solution is correct up to order cfl'^ {p^/"^), since we neglected a term of the 
form q^^'^gQ{q^^'^a). The function go^Q^^'^O') actually is not regular in q: it has an expansion 
in negative powers of q of arbitrary order with coefficients proportional to some power of p. 
This expansion has terms of the form q~^p'' along with terms with higher powers of p and 
same power of g, 

^ l-t^p^qia-"^ 
Fregia,q,p,t) n T^W^- (B.12) 



r{a~^;p 
T{ta-^;p,. 



Here Freg(-) is regular in (/-expansion. If we expand in p to rth order then the term which we 
neglected on the left-hand-side of (B.7) contributes at order ^£ — r. On the other hand the 
second term on the left-hand-side of (B.7) contributes at order + r. Thus, the highest 
order in g-expansion to which the eigen-function are consistent with the assumptions r = ^i. 
We can make several checks of this result. Note that &{ifl) and H can be mapped to each 
other by exchanging t — t- Thus, the (exact) eigenfunctions have to satisfy 

T{ta^^■,p,q)^P{;p,q,t) = ip{;p,q,^) . (B.13) 

This property is obviously satisfied in the above approximation. The approximate eigenfunc- 
tions are orthogonal under the natural measure (up to the order they are valid), 

da T{ta^'^;p,q) 



,P, 



ipeia) ip£'{a) = 5ui ni . (B.14) 



The normalization is given by 



£ = : ni, = l + t, 

i^O : 71^ = 1. (B.15) 

In the special limit p = the approximate eigenfunctions coincide with Macdonald polyno- 
mials up to order i — 1 in q. In the limit p = q = they are identical to Hall-Littlewood 
polynomials (no approximation). Starting from these approximate eigenfunctions one should 
in principle be able to obtain corrections perturbatively in q and p. 



C. Non-maximal punctures 

Let us discuss the procedure of reducing the fiavor symmetry of a puncture by our residue 
calculus. We consider a concrete example that captures the general algorithm. We will derive 
the prescription to partially-close a full puncture to obtain an L-shaped puncture with two 
rows (see figure 7). 
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Figure 7: Closing tlie full puncture down to L-shaped one by closing two minimal punctures. On the 
left we have a generic theory with a superconformal tail ending in an L-shaped puncture represented 
by a Riemann surface. On the right we have a Lagrangian representation of the superconformal tail: 
the squares denote hypermultiplets and circles gauge groups. 



First we couple a theory with a maximal puncture to superconformal tail ending in an 
L-shaped puncture as dipicted in figure 7. The figure on the left corresponds to a Riemann 
surface representation of the theory, and on the right we have a Lagrangian description for the 
tail. The flavor symmetry of the L-shaped puncture is U{1) x SU{N — 2). The superconformal 
tail has two additional factors of U{1) flavor symmetry: in the Riemann surface picture these 
correspond to the two minimal punctures, and in the Lagrangian representation the ^7(l)s are 
flavor symmetries of the hypermultiplets. To obtain the index with the L-shaped puncture 
only, we have to remove the two minimal punctures by computing the residue of the relevant 
U{1) flavor fugacities at t^, as we discussed in the bulk of the paper. We will implement this 
procedure in a convenient way to obtain an explicit expression for this index. The index of 
this theory is given by 



n 

1=1 



dci 
27rict 

n 



Xv(c)A(c)Xc(c-\..) 



(C.l) 



dbi 



2-Kb. 



{b,2/}). 



Xq denotes the index of the Riemann surface to which we glue the superconformal tail. Here 
the free hypermultiplet indices are given by 

x,,p(b,d,z)= n n n 



1 



1 



(djbiz) 



1 



p''+^q"'+^r2djbiz 



i=l j=l m,n>0 
N 

^hyp{c,x,{h,y}) = Y{ n 

i,j=l m,n>0 



1 



p^q^t^djbiZ 



l-pnqnH2[djbiz) 



l-p^+'^q-^+H~2{xbiCj) 



1 



, pn+l^m+l^- 2 



xbiCj 



p'^q'^t^xbiC 



It. J 



pnqm-f-2 (^xbiCj)' 



(C.2) 
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We have defined {61, . . . 1)^} = {ybi^ , ■ ■ ■ ybj/_i,y^~'^}- Let us look for poles which can pinch 
the integration contours when a combination of the U{1) fugacities associated to the U{1) 
punctures is equal to t^. There are three U{1) fugacities in the tail denoted above by x, 
y and z. Let us denote by a and P the combinations of U{1) fugacities corresponding to 
the minimal punctures and by 7 the U{1) fugacity inside the L-shaped puncture. Then we 
have the following identification of the C/(l) fugacities in the two descriptions: the Riemann 
surface and the Lagrangian, 

1 N-2 9 A 

a = x, f3 = -z N ^ ^ = y^zN . (C.3) 

y 

We will consider the poles sX a = fi = one after the other. Note that if the c fugacity is 
not gauged there are no such poles. The poles in the c integration occur at 

1 l-iV 

Ci = — ^<^0) c,- = i = l...iV-l; cn = ^^yy-N ■ (^-4) 

Thus if we choose 

x = t^ , (C.5) 



all the c contours are pinched and we obtain a pole. We evaluate then the integrand of the c 
contour integral at 



We obtain 



x = tl- Cj = y-\^j), j = l...N-l; cn = y^'^ ■ (C.6) 



.N-2 „ 

^7 n ^,Mh-')A{h-')lc{h,..) X (C.7) 
j=i 3 

(b) A(b)X;,,p(b, d, z) T'^yp{{h, y}) . 
The free hypermultiplet which coupled to c becomes 

Iuypic,x,{h,y})^l,^^i{h,y})= [[ [[ . . l-M./l 

i,j=lm,n>0 ^ PI ''Oi/Oj 1 p q Oj/Oi 

= ^ . (C.8) 

Ar!Xv(b-i)A(b-i)Xi/ 
In turn, the index becomes 

, N-2 

2Tib 



n ^Xc{h,..)Iv{\y)A{h)Inyp{h,d,z) . 
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The poles in the b integration contour now occur at 

1 2-jV 

6, = |C'0)' J = l...iV-2; bN-i = l^. (C.9) 
We have 7 = y zn and /3 = jv . Thus we set /3 = 

2-N 

2 



y 

1 ^r_-l _1 N-2 _-, t 2 1 JV-2 

yz = t^-f, 2/^^ 2^ 2 , 67V_iy = 2-/V = ^ ■ (CIO) 



In particular 

1 -I 1 JV-2 1 N-2 

Cj = -dj\ j = l...N-2, CN-i = t2j^, CN = t-^^^. (C.U) 

If we only consider the poles in the decoupled piece (the sphere with the superconformal 
tail) there are no contour pinchings. However, since we know that these poles exist, say by 
decoupling each of the minimal punctures together with a maximal one as before, we deduce 
that Xq has to have a relevant pole. We note that 

TyA(b)Xy(b)X,,,p(b,d,z) ^ / -. (C.12) 

r(t2(7^v/^dj)±^;p,g) 

Thus we deduce that the index (times the two free U{1) factors) is given by 

= xa{ th"^ , rh'^], . . .) I . (C.13) 

7 r(t2(7JV/2^^.)±l;p,g) 

To obtain the index of a theory with puncture corresponding to L-shaped Young diagram 
we have to compute a certain residue of the theory with full puncture and strip off a free 
hypermultiplet. This procedure can be generalized for arbitrary non-maximal punctures. 

Let us compare the above result with the general prescription to compute the index of 
non maximal punctures in the Macdonald case suggested in [23]. Following the prescription 
of [23] the index of the theory with L-shaped puncture is given by 

Xl = Kl{7, d) k{-) Vi-d^^t'^j'^ , r^^) ^a(-) . (C.14) 

A ^ 

Where ^a(") and K{-) are combination of functions not depending on the flavor fugacities of 
the L shaped puncture. The function il)^ is a Macdonald polynomial and we also have^^ 



KLi7,d) = PE 



^ i,j=l ^ i=l 



. (C.15) 



■^^The pletliystic exponential is defined as PE[f{x)] = exp(5]]^-^ jf{x^)) where we also have to specify what 
are the parameters x with respect to which the plethystics is done: in index applications the parameters x are 
all the fugacities appearing in index definition. 
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On the other hand the index of the same theory but with the L-shaped puncture traded with 
the maximal one is given by 



X = K(a) k{-) ^\ai, • • • , ajv) ^x{-) • 



(C.16) 



The factor fC for the maximal puncture has the following form 

N 



)C = PE 



n (t n. I, 



{tai/aj; q) 



(C.17) 



We now want to obtain 1^ from a residue of I. X has poles whenever 

ai = tq^ aj . 

Since Hi^i = 1 a class of poles is 



1 



T-rAT-l 1/2 ■ 



(C.18) 



(C.19) 



Lets us consider the pole at ai with £ = and define bi = Ylf^ti a - Then the fugacity 
assignment in the orthogonal functions will be 



(C.20) 



and of course by construction bf Yl'h^^ ai = 1. The above assignment corresponds to two row 
Young diagram with one box in the first row and — 1 boxes in the second and defining a 
puncture with fiavor symmetry S{U{1) x U{N — 2)). We see that it is the same as in (C.14) 
with the following identification of fugacities, 



JV-2 
6l = 7 2 



ai = dij , 

Let us now compute the residue of the index at this pole 



(C.21) 



1 - 2 1 ^"^ 1 

a,M, ^^""^^ = JMWioKq^) {tai/aj;q) H (t3/2(^./t^)±i. ^)(ii/2(5^/„^)±i. ^) 



(C.22) 



On the other hand we also have from (C.15) 



^^^""^ {t;qW;q) {taja,;q) H (^3/2(^^/5^)^1., 



(C.23) 
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The ratio of the two quantities above is simply given by 



ai 



1=2 



2PE 



^ 

tl/2 / 1 



(C.24) 



Thus the prescription to obtain the index of a theory with S{U (1) x U{N — 2)) puncture from 
the index of the theory with the maximal puncture is just to consider a certain pole of the 
latter discussed here and multiply the residue by the index of a free U{1) vector multiplet 
and divide by the index of an apropriate free hyper- multiplet. We get complete agreement 
with (C.13). The procedure of [23] to compute the index of theories with any non-maximal 
punctures can be phrased as a computation of residues as was discussed in [43] . 
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